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Abstract. Classical (maximal) superintegrable systems in n dimensions are Hamiltonian 
systems with 2n — 1 independent constants of the motion, globally defined, the maximum 
number possible. They are very special because they can be solved algebraically. In this 
paper we show explicitly, mostly through examples of 2nd order superintegrable systems 
in 2 dimensions, how the trajectories can be determined in detail using rather elementary 
algebraic, geometric and analytic methods applied to the closed quadratic algebra of sym¬ 
metries of the system, without resorting to separation of variables techniques or trying to 
integrate Hamilton’s equations. We treat a family of 2nd order degenerate systems: os¬ 
cillator analogies on Darboux, nonzero constant curvature, and flat spaces, related to one 
another via contractions, and obeying Kepler’s laws. Then we treat two 2nd order nondege¬ 
nerate systems, an analogy of a caged Coulomb problem on the 2-sphere and its contraction 
to a Euclidean space caged Coulomb problem. In all cases the symmetry algebra structure 
provides detailed information about the trajectories, some of which are rather complicated. 
An interesting example is the occurrence of “metronome orbits”, trajectories confined to an 
arc rather than a loop, which are indicated clearly from the structure equations but might be 
overlooked using more traditional methods. We also treat the Post-Winternitz system, an 
example of a classical 4th order superintegrable system that cannot be solved using separa¬ 
tion of variables. Finally we treat a superintegrable system, related to the addition theorem 
for elliptic functions, whose constants of the motion are only rational in the momenta. It is 
a system of special interest because its constants of the motion generate a closed polynomial 
algebra. This paper contains many new results but we have tried to present most of the 
materials in a fashion that is easily accessible to nonexperts, in order to provide entree to 
superintegrablity theory. 
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1 Introduction 

Classical and quantum mechanical Hamiltonian systems that can be solved explicitly, both 
algebraically and analytically, and with adjustable parameters, are relatively rare and highly 
prized. Eamous classical examples are the anharmonic oscillator (Lissajous patterns) and Kepler 
systems (planetary orbits) and the Hohmann transfer for orbital navigation [7]. Eamous quantum 
examples are the Coulomb system (energy levels of the hydrogen atom, leading to the periodic 
table of the elements), and the quantum isotropic oscillator. The solvability of these systems 
is related to their symmetry, not necessarily group symmetry. This higher order symmetry is 
captured by the concept of superintegrability. A natural classical Hamiltonian system (with the 
Hamiltonian as kinetic energy plus potential energy) on an n-dimensional Riemannian space is 
said to be (maximally) superintegrable if it admits the maximally possible 2n — 1 functionally 
independent constants of the motion globally defined (usually required to be polynomial or at 
least rational in the momenta). Similarly a quantum Hamiltonian system H = A„+H (where A„ 
is a Laplace-Beltrami operator on an n-dimensional Riemannian manifold and H is a potential 
function) is (maximally) superintegrable if it admits 2n — 1 algebraically independent partial 
differential operators commuting with H. There is a rapidly growing literature concerning these 
systems, e.g., [2, 3, 4, 6, 8, 9, 10, 11, 13, 20, 28, 30, 31, 32, 33, 38, 39, 40, 42]. In this paper we 
consider only classical systems and n = 2 so all of our systems admit 3 independent constants 
of the motion. By taking Poisson brackets of the classical constants of the motion we generate 
a symmetry algebra, not necessarily a Lie algebra, which is never abelian. (This contrasts with 
integrable Hamiltonian systems which admit n constants of the motion in involution, so that the 
symmetry algebra is always abelian.) It is this nonabelian symmetry algebra and its structure 
that is responsible for the solvability of superintegrable systems. 

To be more explicit, along a specific classical trajectory each constant of the motion Cj takes 
a fixed value ij, so the trajectory can be characterized as the intersection of 2n — 1 hypersurfaces 
Cj = ij in the 2n-dimensional phase space. Thus in principle the path of the trajectory can be 
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determined algebraically, though not how it is traced out in time. Since the Poisson bracket of 
two constants of the motion is again a constant of the motion, the nonabelian symmetry algebra 
gives us relationships between the symmetries. In the quantum case each bound state eigenspace 
of the Hamiltonian is invariant under the action of the symmetry algebra, so that a knowledge 
of the irreducible representations of the symmetry algebra gives useful information about the 
dimensions of the eigenspaces and of the eigenvalues themselves. 

In this paper we illustrate the value of superintegrabilty by studying and solving several 
families of classical superintegrable systems. Second degree superintegrable systems are those 
whose generating symmetries are all polynomials in the momenta of degree < 2. All of these 
systems have been classified for 2 dimensions (as have the systems with nondegenerate potentials 
in 3 dimensions) [5, 19, 23]. They occur on constant curvature spaces (with 3 Killing vectors, 
admitting the most symmetries), on the 4 Darboux spaces (admitting 1 Killing vector) and 
6 families of Koenigs spaces (admitting no Killing vectors) [25, 33]. Thus, after the constant 
curvature spaces, the Darboux spaces admit the most symmetries. The Darboux metrics an be 
written as 

D\ : ds^ = 4:x[dx^ + D2: ds^ = ^ ^ {dx^ + dy ^), 

D3 : ds'^ = (dx‘^ + dy‘^), DA{b): ds^ = ^ ^ idx"^ + dy '^), 

sin 2x 

An example of a Koenigs space is 



We first treat some analogs of the harmonic oscillator. (A treatment of the Kepler system and 
its analog on the 2-sphere from this point of view is contained in Chapter 3 of the article [33].) 
We start by studying a system on the Darboux space DA{b) with a 1-parameter potential. 
Then by taking limits (contractions [15, 24]) of this system we obtain systems on the Darboux 
space D3, on the Poincare upper half plane, the 2-sphere (the Higgs oscillator [14]), and finally 
the isotropic oscillator on Euclidean space. 

The second class of examples are nondegenerate (3-parameter potential) systems. The first, 
57 in our listing [23], can be regarded as a caged version of an analog to the Coulomb potential 
on the 2-sphere. It contracts to the caged Coulomb system EIQ in Euclidean space. 

The preceding examples can also be studied analytically via separation of variables. However, 
the Post-Winternitz system [36] cannot, see also [29]. It is 4th degree superintegrable but non- 
separable. However, we show that the classical orbits can be found exactly via superintegrabilty. 

The last examples are different. They are separable in elliptic coordinates and can be derived 
via an action-angle construction. The usual action-angle construction of a superintegrable and 
separable system requires the addition theorem for trigonometric or hyperbolic functions and 
leads to polynomial superintegrability, e.g., [20]. The construction here uses the addition formula 
for elliptic functions [35]. It leads to classical systems where some of the constants are nontrivially 
rational in the canonical momenta. We present two examples where, however, the systems have 
polynomial symmetry algebras, so we consider them as superintegrable and worthy of study. 
The new feature here is that the symmetry algebra closes polynomially, even though the system 
is only rationally superintegrable. These are the first examples of such behavior known to us. 
Again we can determine the trajectories exactly. 

With the exception of the elliptic systems, all of the superintegrable systems in this paper 
have been derived and classified before; we take the structure equations as given and show 
how superintegrability alone leads to formulas for the trajectories. These formulas and their 
analysis are new. The elliptic systems have not been found elsewhere to our knowledge, so we 
demonstrate the procedure to derive them. 
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Figure 1. A diagram showing the contractions of oscillators. 


This paper is partly pedagogical (the Introduction and the Appendix) but mostly new re¬ 
search (Sections 2-5). In the Appendix we give a brief review of fundamental definitions and 
results from classical Hamiltonian mechanics, adapted from [33], needed as background for our 
computations. The point here is to illustrate how, using the structure equations of super- 
integrable systems alone, we can derive and classify the trajectories via algebra. Hamilton’s 
equations are used only to determine the periods of orbits for systems with degenerate poten¬ 
tial. Separation of variable techniques are not used, except in the last section; our approach 
is easy to understand geometrically. Analogously the spectra of the corresponding quantum 
systems can be obtained algebraically from the structure relations, though we do not treat this 
here. Some recent papers, e.g., [26, 27], adopt a related but different approach by using ladder 
operators constructed from the structure algebra to compute trajectories and spectra for sys¬ 
tems with degenerate potential. We know of no prior treatments of the nondegenerate caged 
systems S7 and E16, or of the use of the structure algebra to call attention to special orbits, such 
as those for which TZ^ = 0. We point out the contractions that relate our various superintegrable 
systems. 


2 Examples of 2D 2nd degree degenerate systems 

Our first examples are degenerate systems. These have 1-parameter potentials and always admit 
a symmetry that is a 1st degree polynomial in the momenta, hence a group symmetry that can 
be interpreted as invariance with respect to rotation or translation corresponds to a constant of 
the motion which leads to an analog of Kepler’s 2nd Law. The only possibilities in 2 dimensions 
are constant curvature and Darboux spaces [22]. In [33] there is an example of an analog of the 
Kepler problem on the 2-sphere that satisfies Kepler’s three laws and then by a limiting process 
(a contraction) goes to the Kepler system in Euclidean space. Here we will treat an analog of 
the harmonic oscillator on the Darboux space D4{b) and show that it obeys analogs of Kepler’s 
laws. Then by taking contractions to superintegrable systems on the Darboux space D3, on 
the Poincare upper half plane (equivalent to a system on a hyperboloid in Euclidean space) 
and, finally, to the isotropic oscillator system in Euclidean space, we will see that using ideas 
from superintegrability theory alone we can understand the basic properties of these systems: 
conservations laws, explicit trajectories, etc., and how they are related. Eig. I describes the 
contraction relationships that we will exploit. An important feature of degenerate systems is 
that they always admit 4 linearly independent symmetries that are 2nd degree in the momenta, 
whereas only 3 can be functionally independent. Thus there must be a relation between these 
symmetries. This relation emerges from the structure algebra obeyed by the symmetries. 
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2.1 The D4{b) oscillator 


We consider the superintegrable system 


^ ^ ^ , _a_ 

2cosh(2x) + 6^ 2cosh(2x) + 6’ 

with a basis of constants of the motion, T-L, J = Py, and 

n = - cos{2!,) (pj+pj) _ co.h(2x)p;) 

acos(2y) 

— sinlzT/j smh(2x -^-, 

Vi/; V ^I'xyy 4 (.03^2 _ 2 + 5’ 

( 2 rh(2lT 1 ~ 

asin(2y) 

+ cos(2y) smh(2x)p^py - -- 2 ————, 

4 cosh (x) — 2 + 0 


( 2 . 1 ) 


( 2 . 2 ) 


(2.3) 


Here the spaces D4 are indexed by a parameter b > —2. (In the limit as 6 —)• —2 the space 
becomes a 2-sphere and this system becomes the Higgs oscillator [14].) The variable y can 
be interpreted as an angle, and the space and potential are periodic in y with period vr. The 
constants of the motion generate an algebra under the Poisson bracket obeying the structure 
equations 


{J,'H] = {yi,'H] = {y2,'H] = ^, (2.4) 

{j,Ti} = -2T2, {j,T2} = 2 Ti, {yi,y2} = ^j^+ 2bjn-aj, (2.5) 


with Casimir 


yl + yl = + 'H^+ bJ^'H-aJ^. ( 2 . 6 ) 

Note that there are 4 linearly independent symmetries T-L, Ti, T 2 as polynomials in the 
momenta, but there can be only 3 functionally independent generators. The dependence relation 
is given by (2.6). 

From the first two equations (2.5) we see that (Ti,T 2 ) transforms like a 2-vector under 
rotations about the 3-axis. The sum Ti + TI is expressed in terms of constants of motion so 
the sum is constant. We set Ti + Tf = where k is the length of the 2-vector. Thus we can 
choose a preferred coordinate system such that Ti = « and T 2 = 0. 

To determine a trajectory, we need to express x and y in terms of the constants of the motion 
along the trajectory. We do this by eliminating the momenta from equations (2.1)-(2.6). We 
see that Ti and T 2 can be expressed in an alternate form: 

Ti = — cos{2y)'H -|- cos(2y) cosh(2x)p^ — sin(2?/) smh.{2x)pxPy, (2.7) 

T 2 = — sm{2y)'H + sin(2?/) cosh(2x)py -|- cos(2?/) sm}i{2x)pxPy (2.8) 

To eliminate px, we multiply equation (2.7) by cos(2y), equation (2.8) by sin(2?/) respectively, 
and add them together, which gives us Ti cos(2?/) -|- T 2 sin(2y) = —Ti + cosh{2x)py. Using the 
facts that Ti = /v, T 2 = 0 and J' = py, we get the orbit equation, 

cosh{2x)J'^ — cos{2y)K = 'H. (2-9) 

Since cosh(2x) is an even function in x, there will be two orbits for positive x and negative x 
that are symmetric with respect to x-axis. Therefore, without loss of generality, we restrict 
ourselves to positive x. 
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Analog of Kepler’s second law of planetary motion. We exploit the fact that there is 
a 1st degree constant of the motion. We consider an orbital motion as taking place in the si — S 2 
plane with polar coordinates si = rcosO, S 2 = rsinO where r = (2cosh(2x) + 6)/(2sinh2x) and 
0 = 2y. In these coordinates the metric is 


ds 


2 1 + 26r2 + Vl6r4 + 46r2 + 1^2, 2^o2 

I = 2- 7 - dr + r dO . 

h 4 - ‘)hr2\2 


If a > 0 the potential is attractive to the origin in the plane; if a < 0 the potential is repulsive 
and the trajectories are unbounded. 


Theorem 2.1. Trajectories of the D4{b) oscillator sweep out equal areas in equal times with 
respect to the origin in the si — S 2 plane. 


Proof. In the interval from some initial time 0 to time t the area swept out by the segment of 
the straight line connecting the origin and the object is A{t) = ^ r‘^{9)d6. Thus the rate at 

which the area is swept out is 


dA dAdO 1 2iqi.\\^^ 2dy 


( 2 . 10 ) 


since dy/dt = {dy/dO) ■ (dO/dt) = (1/2) • {dO/dt). To get {dy/dt), we take the Poisson bracket 
of Ti and y, which is 


^ _dy _ 2pysmh.‘^{2x) 
^ dt 2cosh(2x) + 6 


2J"sinh^(2x) 
2cosh(2x) + 6’ 


then plug in this expression for (dy/dt) into equation (2.10), to get, (dA/dt) = J/2 which is 
a constant. ■ 


Analog of Kepler’s third law of planetary motion. 
Theorem 2.2. For a > 0 the period T of an orbit is 

T = \^( 2 + 6 ^ 2-6 \ 

\^(-2U-bUFa) ,/(2n - bFL + a)) ' 


( 2 . 11 ) 


Proof. The total area swept out as the trajectory goes through one period is 



_ ^ ^2, (2 + 6) 2J^'H + J^ + {2- b)y/-K? +'H‘^ - 2jm 
~ 16^ "" V(-k2 + _ 2J‘iU + + 7^2 + 2J^'H + J^) 

However, from the second law we see that A(T) = (T/2) • J and, using the fact that = 
+ Td? + bJ'^'H — aj‘^, the period can be expressed in the form (2.12). ■ 


Now we begin an analysis of the trajectories. We first assume that a > 0 so that the potential 
is attractive. 

Restriction on %. Recall there is a restriction when we try to express r in terms of 0 and 
other constants of the motion. That is, cosh(2a:) = (FL + cos9k)/(J'^) should always be larger 
than 1 for any 9. Then since k > 0, we have the restriction FL — k > implying that for 
a closed trajectory, Fi should always be positive. Here k is the length of the 2-vector (Ti,T 2 )- 
To be explicit, nf = + Fi"^ + bJ^FL — aj^. Then by squaring both sides of an alternate 















Examples of Complete Solvability of 2D Classical Superintegrable Systems 


7 


S,vsS, 



Figure 2. Orbit plot for J' = 10, H = 130, k = 20, & = 0 (equal axes). 


form of the restriction equation, T-L — > k, and plugging in the expression for we get 

(a — (6 + 2)%)j‘^ > 0. Noticing the fact that b + 2, Ti, and j‘^ are all larger than 0, we reach 
the conclusion, T-l < a/{b + 2) for a bounded orbit. For larger values of Ti the trajectory is 
unbounded. 

Case of bounded trajectory. To plot the trajectories on the si — S 2 plane it is convenient 
to write si and S 2 both in terms of 6, which is 


1 lj‘^{2n + 2cose-K + bJ^) ^ 

= 2\j . 

1 /J2(2«+2cos»./v + 6J2) ^ 

(«+oos«.«)2-J4 


(2.13) 


(Note: it is useful to divide the r part of si and S 2 by 77^^ and introduce new constants %' = 'Kj 
and k' = nj so we can eliminate one constant.) Since in si and S 2 , J always appears in the 
form of 77^, so for every positive 77, there will always be a duplicate case for the correspon¬ 
ding —77, and without loss of generality, we can restrict our discussion to 77 > 0. 

A typical plot of the trajectory on the S 1 -S 2 plane is Fig. 2, and as one or more of 7^, J and b 
becomes larger and larger, they will dominate the r term and make r less susceptible to the 
change of cos 0. Thus the plots will look more and more circular. Furthermore, if k becomes very 
large, the plots will tend to move towards the negative si direction and appear in an elongated 
form as in Fig. 3. The plot will always be symmetric with respect to si-axis, because si(0) is 
even and S2(0) is odd in 0. Also, since the plot will always lean towards the negative si side, 
larger r always occurs at smaller si value. 

To trace out an orbit in time, we would need a starting point to integrate Hamilton’s equa¬ 
tions. Conventionally, we choose a point that is closest to the origin, which we call perihelion. 
From the plots of the orbits and an easy analysis of equations (2.13), it is obvious that this 
point is the intersection of the trajectory with the positive si axis. The point on the x-y plane 
that corresponds to this perihelion is (xo,0) where xq = ^ arcCosh('^^). Plugging y = 0 into 
equation (2.8) gives PxPy = 0. Realizing py = J ^ 0, we know px must be equal to 0. We 
see that the perihelion points on the phase space trajectory are uniquely determined by the 
constants of the motion as follows: 

H + K 


cosh(2x) 


y = 0, 


Px = 0 


Py = 77 , 


(2.14) 
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Figure 3. Orbit plot for Y = 1, H = 1002, k = 1000, b = 0 (axes are not set equal). 



Figure 4. S 2 vs si for Y = 10, "H = 120, k = 20, b = 0. 


and we can see that if we know the perihelion point on the phase space as in equation (2.14), 
a set of constants of the motion can be uniquely determined. 

Case of escape velocity. This is the case when Ti = + k, a bifurcation point on 

the momentum map. In polar coordinates, as 0 —)• vr we have r —?• oo, whereas for si, S 2 in 
equations (2.13), si —)• —oo and S 2 —)• 0. A plot for S 2 vs. si is given in Fig. 4. 

In Fig. 4, the two tails will extend to infinity in the direction of negative si-axis, and they 
will get closer and closer to the si-axis for smaller and smaller si value but they will never touch 
the Si-axis. 

Case of unbounded trajectory. Here, there is more than one value of y, i.e., 9/2 that has 
no corresponding real value of x in equation (2.9). 

Example 2.3 {71 < J7"^). In this case, since Ti < , r will blow up before si goes becomes 

negative, so the entire trajectory will be bounded on the positive-si side of the plane. A boundary 
case {% = j"^) is plotted in Fig. 5. For Fig. 5 the constants are J = 10, Ti = 100, k = 20, 
6 = 0, so for 0 < 9 < 7r/2 and 37r/2 < 9 < 2tt, corresponding real values of x exist and there is 
a trajectory. For 9 = 7r/2 and 37r/2, r will goes to infinity and at the same time, S 2 will go to 
infinity and si to 0 which is represented by the two tails in Fig. 5. 

Example 2.4 (J7^ <%< + k). In this case, the trajectory will extend to the negative-si 

side of the plane. A typical plot for an orbit is given in Fig. 6. 
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Figure 5. S 2 vs Si for J' = 10, H = 100, k = 20, b = 0. 



Figure 6. S 2 vs si for J = 10, H = 110, k = 20, b = 0. 


2.2 Contraction to the Poincare npper half plane 

Using the Hamiltonian Ti, (2.1) modified by a constant asH' = T-L — [a/(2 + 6)], we let x = eY, 
y = eX,a = -[{2 + b)P]/[2e% J = /C/e, 3^i - {2J^)/{e^) +%- 2Xi, 3^2 « (2df2)/(e), and go 
to the limit as e —)> 0. Then 


^ 2 (Px + Py) + 


K,= pxi 


- 2XYpxPY - ^ (6 + 2), X 2 = Xp\+ YpxPY + 


The structure equations become 


/3 


^ T 2 , 


{/C,Ti} =- 2 T 2 , {X,X2} = 1C^ + ^{h + 2), {Xi,X2} =-2XXi +——KW, 


and the Casimir is 

,/?(6 + 2) J^^2^ V ^ + ^ir2-i7' 


X 2 + 


+ X^)Xi - = 0 . 


4 


4 


(2.15) 
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Circle 



Figure 7. Circular trajectory on Poincare upper half plane with /3 = 0, 6 = 0, "H' = 2, /C = 1, A 2 = 1. 

If the potential is turned off this is essentially the Poincare upper half plane model of hyperbolic 
geometry; for 6 = 4 it is exactly that. The space here consists of all real points {X,Y) with 
Y > 0. In this limit there is no longer any periodicity. Using these identities to eliminate the 
momenta we arrive at the equation for the trajectories: 

4/C + X - 2 X 2 ^ +fC^{4IC‘^ + p{b + 2))Y^-n'lC^ib + 2) = 0. (2.16) 

Notice that the orbit equation is even in Y, so the trajectory will always be symmetric with 
respect to the Y-axis. When restricting to the Poincare upper half plane with Y > 0, the 
trajectory will be the upper half of the trajectory traced out by the orbit equation. If /3 > 0, so 
that the potential is attractive to the boundary Y = 0, this describes the portion of the ellipse 
+ {y‘^/B‘^) = 1 in the upper half plane, where 

n'{h + 2) 

4/C2 + ,d(6 + 2)’ 

y = Y. 

If /3 < 0 the potential is repulsive. If however, d/C^ + (6 + 2)/3 > 0 the trajectory (2.16) is 
again a portion of an ellipse. There is a special case that when (3 = 0, we will have 
such that the trajectory will be upper half of a circle. A plot of a trajectory when /3 = 0 is given 
in Fig. 7. And an example of elliptic trajectory is given in Fig. 8. 

If AK? + {b + 2,)fj < 0 the trajectory is the upper sheet y > 0 of the hyperbola (x^/A^) — 
{y'^/B'^) = 1 where 

2_ n’{b + 2) 

4/C2 + (6 + 2)/3’ 

y = Y. 

If AKf + (6 + 2)/3 = 0 equation (2.16) becomes degenerate. From (2.15) we now find 
Xi = —Y‘^K? — 2 XYpy1C and have %' = [4Y2/(6-|-2)]-py > 0. Eliminating py from these two ex¬ 
pressions we find the equation {KfY‘^ + Xi)‘^—l-L'K?{b+2)X‘^ = 0 for the unbounded trajectories. 
These are upper halves of parabolas with orbit equation Y^ = — (Ai/ZC) ± ^/T-Lfb + 2)X, and 
we notice that a certain subset olT-L', Ti, /C and b would correspond to two half-parabolas which 
are symmetric with respect to y-axis. Interestingly, when XxjK < 0, the two half-parabolas will 


A2 = 


A'H!K?{bY2) 


x = X- 


(4/C2 + /3(6 + 2))2’ 
4^2 


4/C2 + /3(6 + 2)’ 


A2 = 


AVfK}{bY2) 


x = X- 


(4/C2 + /I(6 + 2))2’ 
4^2 


4/C2 +/3(6 + 2) ’ 
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Figure 8. Elliptic trajectory on Poincare upper half plane with /3 = 2, 5 = 0, "H' = 2, /C = 1, A 2 = 1. 

Hyperbola 


> 3 - 
2 - 

1 


-10 -8 


Origin 


Figure 9. Hyperbolic trajectory on Poincare upper half plane with /3 = —3, h = Ti! = 2, K, = \, 
X 2 = 1 . 


A Parabolic Trajectory 


Axil 


Or^in 


Figure 10. Parabolic trajectory on Poincare upper half plane with (3 = —2, 6 = 0, "H' = 2, /C = 1, 
A-i = 1. 


^ A Intersecting Parabolic Trajectory 

6 - 

5 - .- 

4 - , 

> 

3 - „ 

2 - " 

1 - 

qI-1-1- lA -1-1- 

-15 -10 -5 Origin 5 10 15 

X 

Figure 11. Parabolic trajectory on Poincare upper half plane with /3 = —2, b = 0, H' = 2, 1C = —1, 
Xi = l. 

cross as in Fig. 11. A special case occurs when = 0 so that the trajectory becomes a straight 
line. Since %' = (4y^ • p\)/{b + 2), we see this implies that py = 0 so the trajectory must be 
parallel to the x-axis. Furthermore, if = 0 and Xi/K, = 0, the trajectory will be exactly the 
y-axis. Moreover, for = 0 and Xi/K, < 0, there will be no trajectory. 
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We mention that this Poincare upper half plane model can be mapped to the unit disk, the 
Poincare disk, as well as to the upper sheet of the 2-sheet hyperboloid is 3-space. Further, this 
model can be contracted to the flat space system H = pi + + ocx, E5 in our listing [23]. 

Indeed, if we set 


X = y, 


Y = X + 


Vb 2 

2 e 


/3 


q;\/F+^ 

e 


n{e) 



a(6 +2)3/2 \ 
4e J’ 


then lim ?^(e) = pi +pi + ax. We don’t go into detail here but we give another detailed example 
£—>■0 ^ 
in the next section. 


2.3 Contraction to the D3 oscillator 

The oscillator in the space Darboux 3 (U 3 ) has Hamiltonian 


n 


1 

2 6=^ + 1 


{pI+pI) + 


13 

+ 1 


(2.17) 


with a single Killing vector J = Py and a symmetry algebra basis, {?/, Yi, T 2 }. Here, 


^2x 


^ 1 X ■ 1 e- 

Yi = -e sinyp^py + 


cos yp^ 


^2x 


V lx ,1 e 

/f2 = --e cosyp,py + --^^ 


1 J-2) 2 1/3cosy 

4 + 1 2 e^TT’ 

2 1 -h 2) . 2 , l/3siny 

smyp^ - - ^ ^ ^ smypy J- -- 


4 -k 1 


2e’’^-kl 


The structure relations are 

{j7,Yi} = -Y2, {j7,Y2} = Yi, {A:^^x2} = ^jn-^j, 

and there is the functional relation + X 2 — J- § = 0. 

Theorem 2.5. The D3 oscillator system is a contraction of the DA oscillator system. 

Proof. We can get (2.17) as a limiting case of (2.1) by taking 

x = |--^ln(e), = 31' = Sen, “ = (^-1^) 

Then we have %' = 5 Y p\i), J' = 2J = pyi, and y[ Ri -4eTi, T 2 ~ - 4 eT 2 - 

First, it is obvious that J' = 2J = pyi. Then as e —)■ 0, since 

g 3 :'-ln(e) _j_ ^-x'+ln(e) gx' _j_ g 2 g-x' gx' _ g 2 g-x' 

cosh(2x) =- - -= - — -, sinh(2x) =- — -, 

plugging these two identities into the Hamiltonian equation (2.1) and also using equation (2.18), 
we have 


2,^4,^-2x' 


n' = 8 e 


-2e'^-|-e'^e 


“iW 




+ 




h 


+6^e ^ I 1 
e ' e ‘ 


1 e^^' - 2e2 -k 


2 -k en-^ -k 


— (P^/ YP^yl) + 


/3 


1 e 


+ e2e-’"' + 1 2 e’"' + 


2x' 


+ 


13 

e^ + Y 
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because the terms with e in front of them will all go to 0 as e —)• 0. This %' is exactly the 
Hamiltonian of the Darboux-3 oscillator. 

Similarly, by substituting the identities in equation (2.18) into equations (2.7) and (2.8), we 
have. 


r, = -H- COS!/ + (jP-' + -iPy) - "“P 


/3 


+ COS y 


e^' + e^e 

2 ^ 


1 


X 

2 ^ 


8e(e^' + e^e + 1) 

1 1 

iP-'Py' = 


1 / 2e^ (e^ + 1) e 

4 I ^ 


2 x' 


1 


e^' + 1 


/ 2 ,1 x' ■ / /^ 1/3cosy' 

cos(y)py+ -e s\n[y)p^,py, + --^r^ 


1 e 


2x' 


, 2 le"=(e"^ + 2 ) ^ 2 , 1 x' w , 1/3cosy' 

- 4 gx' + 1 )Px' - 4 gx' + 1 ^<y ) Py' + 2"" ^P-'Py' + 

as e —>■ 0. In the same sense, 

^,/ 1 2 le"^'(e"^'+ 2) 2 1 x' / , l/3siny' 

^2 = 4 ^ 7 ^ sm(y )p,, - sm(y ) p^, - -e cos(y )p,/yj,/ + 2 ^ 7 ^- 

We notice they are exactly the same vectors as Xi and T’ 2 . So we have proved that in this 
limiting case, Darboux-4 space contracts to Darboux-3. ■ 


Using this transformation, we can get the orbit equation for the D3 oscillator from the orbit 
equation for the D4(6) oscillator. First, since k'"^ = y'\ + y'\ = 166^(3^1 + ^ 2 ) = Ibe^K^ and 
assuming e > 0, we have k’ = 4eK. Then putting k = and the identities (2.18) into (2.9), 


we derive 


+e^e 


2e 


■) (" 2 “)^ “ cos(y')(|^) = Thus, e^' - 2cos{y')n' = 77' as e —> 0. 


2.3.1 Embedding 133 in 3D Minkowski space 


We can embed D3 as a surface in 3D Minkowski space with coordinates X, Y, Z in such a way 
as to preserve rotational symmetry. For example, let 


X = 2\/2e"tyr 


■ y 

cos 


Y = 2V2e-^ + e-- 


: . y 

sm 


Z=^ln 

12 


/#(6 + 5 e"^)V 3 + 2e2^ + 5 e^ + l\ ^- 

-p-^- - e""^V2V3 + 2e2^ + 5e^. 

V^(6 + 5e^)V3 + 2e2^ + 5e^ - 1 / 


Then dX'^ + dY'^ — dZ^ = ^^ 43 ^(dx^ + dy^). Such embeddings are not unique. 

Discussion of trajectories. To see how the trajectory behaves on this surface, we impose 
the orbit equation e^' — 2cos(y')K' = 77' to X, Y and Z. From the orbit equation, we have 


77 + 2 cos y • k 
772 


7P 


+ 


2 k 

7p 


cosy = a + 6 cosy, 


(2.19) 


where we ignore all the primes on the letters for convenience. Also, we treat a = ^ and b = ^ 
as two new constants, so the number of constants is reduced by 1. Then we can express X, Y 
and Z all in terms of y. We notice that since > 0, the trajectory will be closed when a > b 
and unbounded when a < b. 


Example 2.6. Plots of an “elliptical” shape trajectory and its overhead 2-D view are given in 
Figs. 12 and 13. (The size of the “cone” that appears on the graphics is determined by the 
minimum value we choose for e® in MATLAB. The smaller the minimum value, the larger the 
“cone”. For these two plots, the minimum value for e* we choose is 1 and we adjust to different 
sizes of “cones” for other examples to make the plots clearer.) 
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D3 EtnCiedillng 



Figure 12. 3-D View of trajectory with a = 2, 6 = 1. 



Figure 13. Overhead view of trajectory with a = 2, 6 = 1. 

Example 2.7. As a becomes smaller, the “ellipse” is prolonged, see Figs. 14 and 15. Seen from 
above the shape of the trajectory two merging ellipses. (The minimum value for e* we choose 
is 0.1 for these two plots.) 

Example 2.8 (escape velocity). For a = b, we encounter the boundary case between bounded 
trajectory and unbounded trajectory. A plot of the case is given in Fig. 16. The shape of the 
trajectory when looking from above is like two very “thin” parabolas meeting at the top of the 
“cone”. (The minimum value for e* we choose is 0.001 for these two plots.) 

Example 2.9 (unbounded trajectory). When a < b, there are no real x values corresponding 
to a range of values of y, clearly seen from equation (2.19), so the trajectory is unbounded. 
Examples are given in Figs. 17 and 18 which are a standard 3-D view and a closeup overhead view. 

2.3.2 Contraction of the D3 oscillator to the isotropic oscillator 

The isotropic oscillator in Euclidean space has the Hamiltonian 

= pI + pI + {x^ + 2 /^) 


( 2 . 20 ) 
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Figure 14. 3-D View of trajectory with a = 1.1, 6=1. 



Figure 15. Overhead view of trajectory with a = 1.1, 6 = 1. 


D3 Embedding 



Figure 16. 3-D View of escape velocity trajectory with a = 6 = 1. 


and basis symmetries 

ic = xpy- ypa^, ^ ^ 

The structure equations are 

{£i,/C} = 2£2, {£2,/C} = -2£i, 


{£i,£2} = -2w2/C. 
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X 


Figure 18. Closeup overhead view of trajectory with a = 0.7, 6 = 1. 


The functional relation is C\ + C\ = \T-L^ — uP'K? . We can obtain this system as a limit of the 
D3 oscillator as follows: In (2.17) we set x = x' + In(^), y = y', %' = 2e'}i, j3 = (2a;^)/e^. Then 
as e —>■ 0 we have K = pyi and 


+pI') + 

! 1 // 2 2 w^cosy' 

= eYi = — \ sm.yp^,py, + - cosy [p^, - Py,) 1 H-—, 

V' Y ' ,1 ■ // 2 2 , w^siny' 

^2 = eX 2 = — {-cosyp^'Py! + -smy [p^, - Py,) 1 H -—. 


In terms of flat space Cartesian coordinates X = rcosO, Y = rsin0 we have 

e"' = ^, y' = 2 e, W = pI+pI + u‘^[x^ + y^), 

= ^ (Px - Py) + - y^), ^2 = ^PXPY + ^u;^XY, 










Examples of Complete Solvability of 2D Classical Superintegrable Systems 


17 


a<2 

4 1 -,-,-,-,-,-.-,- 1 

3 . 

2 ■ 

1 . 

>N 0 ■ 

-1 ■ 

-2 - 

-3 - 

-4 I_I I I I I I I_ 

^- 3 - 2-101234 

X 


Figure 19. Elliptical trajectories. 


with /C = 2 {XpY — Ypx ) • From the expressions of %' and X [, we get 


2 W + 4X( - 2uj^X^ 

pI = - 


2 W - 4X( - 2a;2y2 
Py = ^ - 


2 ' 2 
Then from the expression of X 2 , we have PxPy = 2 X 2 — uP'XY. Squaring this equation and 
equating it to the product of the two equations ( 2 . 21 ) we obtain the orbit equation, 

2u? {%' -4X[)- IQuj‘^X' 2 XY = %' - 16Xf - l&X'^. 

The shape of the orbit depends only on SX 2 /{ 1 -L' — 4X(), so we are really investigating the 
equation X‘^ + Y‘^ — aXY = 1 where a = 8 X 3 /{%' — 4X[). We take the right side to be 1 since 
we now only care about the shape, not the scale. 

When 0 = 0, i.e., X 2 = 0, the trajectory is just a circle. When 0 < a < 2, i.e., 0 < 
8 X 2 /(Ti' — 4X[) < 2, the trajectory is a tilted ellipse with y = x as the major axis. In Fig. 19, 
we show three trajectories for different values of a, and also the major axis. As a increases, the 
ellipse becomes more elongated in the major axis direction. 

When 0 = 2, i.e., %' — 4X( = 4 X 3 , there is a bifurcation point on the momentum map. 
The trajectory splits into two parallel straight lines. For o > 2, i.e., 8 X 2 /{T-L' — 4X() > 2, 
the trajectories are hyperbolas with y = x as the symmetry axis. In Fig. 20, we present 
three hyperbolic trajectories, shown in darker colors as a increases. It is also possible that 
8 X 2 /{%' — 4X() < 0, in which case the trajectories will be symmetric about the line y = —x. 


2.4 The Higgs oscillator S3 

The classical system S3 on the 2-sphere is determined by Hamiltonian [14] 


( 2 . 22 ) 


where J'l = S 2 P 3 — S 3 P 2 and J 2 -, Jz are cyclic permutations of this expression. For computational 
convenience we have embedded the 2-sphere in Euclidean 3-space. Thus we can write 

PL + {sipi + S2P2 + saPs)^ 


n// 2 , 2 , 2 , ® 

^ = Pi +P2+P3 + ^ = 


si + si + si 
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a>2 



Figure 20. Hyperbolic trajectories. 


3 

and use the Euclidean space Poisson bracket {J^,G} = + dp.TdsiG) for our 

i=l 

computations, but at the end we restrict to the unit sphere: sf + + '®3 = f sipi + S 2 P 2 + 

■53^3 = 0. The Hamilton equations for the trajectories Sj{t), Pj{t) in phase space are 

^ = {W.»j}. i = 1.2.3. 

The classical basis for the constants of the motion is 
2 

Ti = J'2 + a%, C2 = JiJ2-a^-^, X = Jz. (2.23) 

«3 «3 

The structure relations are 

{T, £ 1 } = -2£2, {T, £ 2 } = 2£i - ?^ + + a, {£ 1 , £ 2 } = -2(£i + a)T, (2.24) 

and the Casimir relation is 


L\^L\- CiU + £ 1^2 + + aCi = 0 . 


To analyze the classical trajectories it is convenient to replace the basis elements £ 1 , £2 in the 
algebra with the new basis set 


5i = 


1 

2 





(2.25) 


Thus 5i = £1 + \{X‘^ — T-i + a), S 2 = £ 2 - Now the first two equations (2.24) become {J73,5i} = 
—252, {J73,52} = 25i, so (5i,52) transforms as a 2-vector with respect to rotations about the 

3-axis. Indeed, a rotation through the angle /3 about the 3-axis rotates the vector by 2(3. The 
remaining relations now become 

{ 5 i,52} = -X{n- X'^ + a), Sf +Si = ^{X^ -n + a)^ - 
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This verifies that the length k of the 2-vector is unchanged under a rotation and shows that 
(5i, 52 ) is similar to the Laplace-Runge-Lenz vector for the Kepler problem in Euclidean space. 
(However, a better superintegrable 2-sphere analog of the Kepler problem is the potential 
ass/ s\ + S 21 called 56 in our listing [23].) In analogy with the choice of periaptic coordi¬ 
nates to simplify the Kepler problem, we choose the preferred coordinate system such that the 
vector (5 i, 52) points along the 1-axis, i.e., 52 = 0, 5i = k > 0. Then equations (2.22) and (2.23) 
become 


J1J2 — 


OiSlS2 


Cf\ — ^2 — T 


a{sl-sf) 


Ji = 


Ji + Ji + Ji 


= n-^. 

4 


(2.26) 


The last 3 equations can be solved to give 


Ji = K + ^{n - - a) - 

Z S 3 

Ji = -Kx\{n-x^Jl = x\ 

Z S 3 


Substituting these results into the square of the first equation (2.26) and simplifying, we get the 
result (assuming a / 0) 


n + 


n-x^- 


a 




-K -I- 


n-X^-a 


s2 - = 0. 


(2.27) 


This is the equation of a cone Asf Bs^ + Cs\ = 0. The trajectories lie on the intersection 
of this cone and the unit sphere sf -|- s| -|- s| = 1. Thus we get conic sections again, as in the 
Euclidean Kepler problem, but this time the sections are intersections with the unit sphere, 
rather than planes. The possible types of trajectory will depend on the signs of A, B, C. 

The projection on the si — S 2 plane is the curve 


n + x^-a 


+ K 




% + X - a ^2 v2 
--- K ] S 2 = X . 


(2.28) 


The identities 


K -b -(77 — X^ — a) 

—K + -(77 — X^ — a) 

= 

(2.29) 

K ~ 

— K 2 ~ 

= 77^2, 

(2.30) 


— K 2 ~ 

= q;77 

(2.31) 


will be important in the analysis of trajectories to follow. 

Analog of Kepler’s second law of planetary motion. We see from equations (2.28) 
and (2.30) that for the nonzero angular momentum X the projection of the motion on the unit 
circle in the si — S 2 plane is a segment of an ellipse, hyperbola or straight line and that none of 
these curves pass through the center (si,S 2 ) = (0,0) of the circle. Thus as the particle moves 
along its trajectory (si(t), S 2 {t), S 3 (t)), the line segment connecting the projection (si(t), S 2 {t)) to 
the center of the circle sweeps out an area. Introducing polar coordinates si{t) = r{(p{t)) cos 
S 2 {t) = r{4>(t)) sin (p(t) we see that in the interval from some initial time 0 to time t the area 
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swept out is A{t) = ^ Thus the rate at which the area is swept out is ^ = 

Now note from Hamilton’s equations that along the trajectory 

X = J‘i = Sl{t)p2{t) - S2{t)pi{t) = Si^ - S2^ = 

Thus ^ Y, a constant. 

Theorem 2.10. If the angular momentum is nonzero, the projection of the trajectory on the 
unit circle in the si — S 2 plane sweeps out equal areas in equal times. 


Now we begin an analysis of the trajectories. There are two cases, depending on whether the 
potential is repulsive (a > 0) or attractive (a < 0). 

Case 1: a > 0. This is the case of a repulsive potential. The equator of the sphere repels 
the particle. Thus motion is confined to a hemisphere. We start an analysis of the types of 
trajectories. 

If Y 7 ^: 0 and using the fact that all trajectories will be periodic for a repulsive force, we see 
that there will necessarily be at least one point on the phase space trajectory for which p 2 = 0. 
Substituting into the phase space conditions (2.26) it is easy to see that this is possible only 
for Si = 0. Solving all the equations completely we find that these points on the phase space 
trajectory are uniquely determined by the constants of the motion as follows: 


P2 = 0, P3 = 0. 


Pi — —n + — {fH, + — cr), 


(2.32) 


This prescription gives us a point on each trajectory, comparable to the aphelion for the Kepler 
system, where we can start to trace out the orbit. It remains to analyze the possible orbits. To 
see what is the available parameter space for the case X 0, note that first we must require 
«: > 0. Noting the identity (2.29) we see that the quantities in brackets must have the same 
sign. However, if that sign is negative then the cone (2.27) will degenerate to a point and not 
intersect the sphere. Thus to obtain trajectories it is necessary that the constants of the motion 
satisfy —k + ^("H — — a) > 0, K > 0. On the other hand, if these conditions are satisfied, we 

see from equations (2.32) that there exist trajectories for which the corresponding constants of 
the motion are assumed. Thus the conditions are necessary and sufficient. 

An Analog of Kepler‘s 3rd law of planetary motion. For nonzero angular momen¬ 
tum X the projection of the motion on the unit circle in the si — S 2 plane is an ellipse (2.28) 
enclosing the center of the circle. The area of this ellipse is easily seen to be ttX / ^/7^. Let T be 
the period of the trajectory. Thus T is the length of time for the projection of the trajectory to 
trace out the complete ellipse and A{T) = ttX /y/lf. Since the area of the ellipse is swept out at 
the constant rate ^ = ’f we have A{T) = Equating these two expressions for A{t) 

and solving for T we find T = 2'k/\/I-L. 

Theorem 2.11. For X ^ 0 the period of an orbit is T = ‘Irx 

If A = 0 then from (2.27) and (2.28) we see that the projection of the motion in the si — S 2 
plane is the straight line segment si = 0. Thus this motion takes place in a plane and the 
trajectory is a portion of a great circle passing through the poles of the sphere. Here, k = 
\{I-L — a), % > a > t). The motion is periodic and the points of closest approach to the equator 
are those such that s‘^ = {FL — 0)1%. The period is again T = There is a special case 

of equilibrium at a pole when FL = a. 

Case 2: a < 0. This is the case of an attractive potential, where the equator of the sphere 
attracts the particle. Again, motion is confined to a hemisphere. We first consider the case 
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where 7 ^ 0. Without loss of generality we can assume Af > 0. Then the projection of the 
motion is counter-clockwise. Now k > 0 and from the identity (2.30) we see that there will be 
3 classes of trajectories, depending on the value of 71. 

If Ti > 0 then the coefficients of sf and in the projection formula (2.28) must have the 
same sign, necessarily positive for a real trajectory. Thus 

H> 0 : n + X'^-a>2K>%-X‘^-a, 2 k+ % - X'^ - a > Q, 

and the projection will be a segment of an ellipse. It is straight-forward to check that the ellipse 
always intersects the circle at 4 points. All of these trajectories lead to annihilation at the 
equator. There will necessarily be a “perihelion” point on each trajectory: S 2 = 0, pi = ps = 0, 
sf = '[k + + X'^ — a)]^, P 2 = K + \{'H + X"^ — a). 

If < 0 then the coefficient of sf is positive and the coefficient of is negative in the 
projection formula (2.28). Thus 

< 0: % + X"^ — a — 2 k < 0, 'H + X^ — a -|- 2 k > 0, 

and the projections will be segments of hyperbolas. Each branch of the hyperbola intersects the 
circle at exactly 2 points. Again, these trajectories lead to annihilation at the equator. Again 
there is a “perihelion” point on each trajectory: S 2 = 0 , = ps = 0 , \^K+\{Ti.+X‘^—a )\, 

P2 = k+]^{'H + X"^ — a). 

li% = Q then k = \{X‘^ — a) so the coefficient of sf is positive and the coefficient of is 0 
in the projection formula (2.28). Thus 

H = 0 : X‘^-a = 2K 

and the projections will be segments of straight lines, sf = X‘^/{X'^ — a). Here pi = 0 and 
P 2 = X‘^—a, a constant in accordance with the analog of Kepler’s Second Law. These trajectories 
lead to annihilation at the equator. There is a “perihelion” point on each trajectory: S 2 = 0, 
Pi=P3 = 0, sf = X^I(X^ - a), pI = X'^ - a. 

Einally, we suppose A = 0, so that the motion takes place in a plane through the poles. 
If — a > 0 then the motion takes place in the plane si = 0. All trajectories annihilate at 
the equator. The projected trajectories each pass through the center of the circle, and we have 
P 2 = % — OL — T-Ls^. If H — a < 0 then the motion takes place in the plane S 2 = 0. All trajectories 
annihilate at the equator. The projected trajectories do not pass through the center of the circle, 
and we have p\ = T-L — a — T-Ls\. li T-L = a then there are possible trajectories on any plane 
passing through the poles. If we go to new rotated si, S 2 coordinates such that the motion takes 
place in the plane si = 0, then pi = 0 and p^ = —as^. All trajectories annihilate at the equator, 
except for unstable equilibria at the poles. 


2.4.1 Contraction of D4{b) to the Higgs oscillator on the sphere 

By letting e —)■ 0 in 6 = —2 -|- and defining new constant a = 4/3, the Hamiltonian equa¬ 
tion (2.1) becomes 

(2sinhxcosha;)^ . 2 2 'i 


n = 


sinh^(2x) / 2 4/3 

ipl+pl) + 


2 cosh(2x) — 2 


2 cosh(2x) — 2 (2 sinh xY 


■ {pi + pD + 


4/3 

(2 sinh x)^ 


- cosh2 x(n2 + , /3(cosh^x- sinh ^ 

-cosh x[p,+py)+ 

We can ignore the constant term —/3 and write 

^. 9/9 9 \ /3 cosh^ X 

n = coslY x{pI + pI) + . , 2 —’ 

smh X 

such that 1-L' = 71 + f3. This is the Higgs oscillator on the 2-sphere. 


1 9 / 9 9 \ cosh^ X 

= cosh X.+P,) + ^r^-l3- 










22 


Y. Chen, E.G. Kalnins, Q. Li and W. Miller Jr. 


In terms of the Euclidean embedding of the sphere we have si =cosy/coshx, S 2 = siny/cosha:, 
S 3 = sinhx/coshx, sosf+s|+s| = 1. If we set si = r' cos 6', S 2 = r'sinff' as in polar coordinates, 
then r' = 1/cosh x and 9' = y. However, we notice that in the analog of Kepler’s 2nd law part 
of the Section 2.1, r = y^2 cosh(2x) + b/ (2 sinh 2x) and 9 = 2y. As e —^ 0, r = 1/(2 cosh x). 
Thus r' = 2r and 9' = 9/2, so this contraction space is a double covering of the original DA{h) 
space. We should be careful of this fact in the further computations. Expressed in the phase 
space (si, S 2 , S3,pi,p2,P3), we have %' = Ji + J 2 + Ji + 4, where Ji = S 2 P 3 - S 3 P 2 and J 2 , J 3 
are cyclic permutations of this expression. Considering the potential part of the Hamiltonian, 
we have a new basis set, y'l = j/ - - fi{s\ - sD/sl, y '2 = 2 J 1 J 2 - (2/3siS2/s|), J' = J 3 . 

We will show that in this limit, the orbit equation (2.9) will yield the orbit equation (2.28). 
Eirst, we transform cosh(2x) and cos(2?/) in equation (2.9) to express them in terms of si and S 2 . 
(s 3 can be expressed in terms of si and S 2 -) From the fact that cosh^ x = l/(sf + s^)) we have 

2 s'^ - 3 “^ 

cosh(2x) = 2 cosh^ x — 1 = ^ — 1, cos(2?/) = cos^ y — sin^ y = -1. 

Si + S2 Si + S2 

Then equation (2.9) becomes [2/(sf + s^) — l]i7'^ — [(si — S 2 )/(sf + s^)]k' = "H, where k' is the 
length of the two vector (T^i,T^ 2 ) and we set y '2 = 0 such that k' = y'l. Multiplying the 
equation by (sf + s|), plugging in the fact that T-L = T-L' — fi and doing some rearrangements, we 
obtain 


W k' 


s? + 


%' yj''^ -p- k' 


si = J'\ 


(2.33) 


where we realize that k' = y'l which is 2 times the basis in equation (2.25), so k' = 2k 
where k is the one in equation (2.28). Therefore, we can see that (2.33) is the same orbit 
equation as (2.28). 

Analog of Kepler’s second law of planetary motion. Following the same procedure 
as in Section 2.1, and noticing that r' = 2r and 9' = 9/2, we have, ^ ^ (4r^(0(t))) ^ = 

= 2^. Therefore, = J', which matches the expression for the Higgs oscillator. 
Period of the orbit. With the same procedure as in Section 2.1 and r' = 2r and 9' = 9/2, 
we find area swept out as a trajectory goes through one period T' is A' = ^ (9')d9' = 

2 Jq'' r'^(9)d9 = 4A where A is the area in equation ( 2 . 12 ). Also, from the second law we see 
that A' = ^J'T'. Hence, 

A' 4A ^ ( 2 + 6 2-6 \ 

~ J' ~ J' ~ ~ y ^{-2n-bnya) yj{2n -bnya)) ’ 

where T is the same as equation (2.11). 

In the limit as e —)> 0, we take 6 = —2 into the above equation of T' and get , 

^ 4 ^ _27r_^ TT 

^mya - /3) + 4/5 Vw' 

which differs from the period in Theorem 2.11 by a factor of due to the fact that the con¬ 
traction is a double covering of the original U4 space. 


2.4.2 Contraction of the Higgs to the isotropic oscillator 

This contraction has a simple geometric interpretation, the contraction of a 2-sphere to a plane. 
We can consider the Higgs oscillator as living in a 2-dimensional bounded “universe” of radius 1 
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in some set of units. Suppose an observer is situated in this universe “near” the attractive north 
pole. Our observer uses a system of units with unit length e where 0 < e <C 1 and we suppose 
that using these units the universe appears flat to the observer. Thus in the observer’s units 
we have si = eX, S 2 = eY, S 3 = + Y'^) = 1 — {X^ + Y'^)e^I2 + O(e^). Here, is 

so small that to the observer it appears that S 3 = 1. Thus, to the observer, it appears that the 
universe is the plane S 3 = 1 with local Cartesian coordinates {X,Y). We compare the actual 
system on the 2 -sphere with the system as it appears to the observer and we assume that is 
so small that it can be neglected, unless we are dividing by it. Now we have 

VX Vv 

Si = eX, S 2 = eY, S 3 « 1, pi = —, p 2 = —, P 3 ~ -{Xpx + Ypy). 

We define new constants h hy a = E — — uPje'^ = h/e^, where J 3 = X. Then 

substituting these results into the Higgs Hamiltonian equation T-L = + J 2 + J 3 X otjs\ = E, 

we find {p\ +Py)/e^ (e^y/X^~+^) = h/e^. Multiplying both sides of this equation by 

we obtain the Hamiltonian equation for the Euclidean space isotropic oscillator 

-H^p\+pl + u:\X^ + Y^)=h, 

in agreement with (2.20). Using the same procedure we find that the constants of the motion 
become fC = X = xpy — ypx, and 

Cl = lim e^£i = Py + £2 = lim e^C 2 = —pxPY — cu^XY, 

an alternate basis for the symmetries of the isotropic oscillator. Similarly, the structure equations 
for the Higgs oscillator go in the limit to the structure equations of the isotropic oscillator. 


3 Examples of 2D 2nd degree nondegenerate systems 


3.1 The system S7 on the sphere 

The classical system S7 on the 2-sphere [23] is determined by the Hamiltonian 


n 


Ji + J2 + Ji + 


aisi 

sW4 + 4 



Q 3 S 3 

V'S? + ^2 


where J'l = S 2 P 3 — S 3 P 2 and J 2 , Jz are cyclic permutations of this expression. We have embedded 
the 2-sphere in Euclidean 3-space, so 


2 I 2 I 2 '7’2 I /t’2 

£1 + £2 + £3 ~ +■ 1^2 




3 

and we can use the Poisson bracket {F, G} = X {~9siFdp.G + dp^Edg^G) for our computations, 

i=l 

but at the end we restrict to the unit sphere: s\ + s\ + s\ = 1 and sipi + S 2 P 2 + S 3 P 3 = 0. The 
Hamilton equations for the trajectories Sj{t), Pj{t) in phase space are 

^ = ^ i = 1.2.3. 

The classical basis for the constants of the motion is 
r _ ^2 aiy/sf + sisi , a2[sl + sl) 

L,i — e/3 2 + „2 ’ 

*2 *2 
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^ _ n- n- , ai^3('S2 + 2Si) 025153 , OS-Sl 

'^2 e/le /3 — 2 A~2-2 ' —~ 


252 V Si + S 2 


2-\/ 5^ + s\ 


The structure relations are TZ = {£i, £ 2 } and 
{TZ, Cl} = 0103 + 4£i£2, 

2 2 

{7^, £ 2 } = -6£? - 2Cl + 4^£i - 2 o 2 ^ + 4o2£i + 

7^2 + 4£? - ACln + 4£i£i - 4a2£? + 4o2^£i - a\% 

— (o| — o^)£i + 2 ai 03£2 + 02 a| = 0. 

3.1.1 S7 in polar coordinates 

In terms of polar coordinates r, 9 where 5i = r cos0, 52 = r sin0, 53 = ±\/l — r^, and 0 < r < 1, 
0 < 0 < TT, we have 


^ ^ 

^ 9 2ai cos 9 + 2 o 2 

Ci=pt + — - 

1 — cos( 2 t^) 


2 X 2 OiCOS0 02 03 \/l — 

jPr ^~9-9 /i~^9 -2 /i^ ’ 

sin sin 0 r 


£2 = ±Vl -r' 


i: f cos9p0 . n \ , /; - rfoi + oi cos^ 0 + 2 o 2 Cos 6 *) 1 

^ ' - h sin 9pr ] pe ± \/ 1 — -'— 


2 r(l — cos^ 9) 


+ -03 cos^. 


Here, ± is interpreted as + in the northern hemisphere and — in the southern hemisphere. 
Hamilton’s equations give 

2n a 2p0 . oi + 2 o 2 cos 0 + Ol cos^ 0 

r = 2(l-r)pr, 6'=^, p0 = - 9 ■ 3fl -> 

sm 9 


2 2oicos 9 202 


Pr — ~3“ + 2rpj. H i ■ 2 a ^ 1 ■ 2 0 ^ / 9 

oTn-^ H T’O cmT^ H * /I _ <r>2 


sin^ 9 
03 \/l — 


Assumptions. We require 03 < 0, and 01,02 > 0. and restrict our attention to trajectories 
caged in one hemisphere (eastern or western): 

— 1 < 53 < 1, 0 < 52 < 1, —1 < 5 i < 1, or 1 > r > 0 , TT > 0 > 0 . 

Note that the north pole is attractive and the south pole is repulsive. 

3.1.2 Trajectories 

Eliminating pg, pr in the expressions for PL, £ 1 , £2 we find the implicit equation for the trajec¬ 
tories: 

0 = cos^ 0o|r^ — 4£^r^ cos^ 9 + A'Hr'^Ci cos^ 9 — 4 cos 9asr‘^C2 — 8 cos 9CiC2r^/l — 

— 4£ioi cos 9r‘^ — 2oi cos 9azriJ 1 — + 4'Hr‘^ai cos 9 — r‘^a\ + 4C\ — 4 £i 02 

— 4oi£2r’\/l — + of + 4£2r^ — 4o203\/l — r^r -|- 4'Hr^a2 — 4'Hr‘^Ci 

-|- 4o3\/1 — r^rCi- 

This is a quadratic equation for cos0 as a function of r, with solutions 


cos 9 = 


A(r) ± 2 ^;^ 
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in the upper hemisphere, where 

N = + 2a3r£,2 — ^T-Lrai + 4/li£2'\/l — + 2rCiai, 

s= {nr^ - ^/l-r^ras- Ci)TZ^, D = r{mCi - ACj + aj). 

We note that 

— 45 = (og + ATiCi — 4£^) (4£2^^ — 4a3-\/l — r‘^ra2 + 4£i03-\/l — r^r + 4'Hr^a2 
— 4'Hr‘^Ci — 4£ia2 — — 4ai£2'\/w-^?’ + 4£i + of). 

In the lower hemisphere the trajectories are given by 


iV(r)±2y^ 

cos d = -=— -, 

D{r) 

N = —aia^^/l — 4- 2a^rC2 — 2'Hrai — ACiC 2 \/l — + 2r£iOi, (3-3) 

s= {nr^ + ^l-r^ras- Ci)n^, D = r{AnCi - + al). 


Since TZ^ > 0 we see that if ^ < £i then 5 < 0. Thus for this case, the trajectory is never in the 
lower hemisphere. If > £i an elementary analysis shows that 5(r) >0 exactly in an interval 
0 < a < r < 1 and is nondecreasing on that interval. 

Case 1: 02 > ai > 0. All trajectories are closed and periodic. We must have £1 > 0 for 
trajectories and we initially assume TZ^ > 0. The equations for perigee and apogee (distance 
from the projection of the trajectory in the equatorial plane to the origin) are 

2 _ a2 + 2nCi ± as^aj + AnCi-ACl 
2(ai+H2) • 

Thus, in order to have physical trajectories we must have 
a| + 4?^£i — 4£g > 0 and a| + 2?^£i > 0. 


Suppose the trajectory touches the equatorial plane at r = 1, 0 = 9q. Then taking a limit in (3.1) 
as the trajectory goes to the boundary we find that if the trajectory touches the equatorial plane 
it does so at an angle cos 0 which is a solution of the quadratic equation 

— £i)£i^ cos^ 9 — {azC 2 — ai{T-L — £1)) cos6* + C 2 — {T-L — £i)(£i — 02) = 0. 


Indeed, cos 6*0 = [2(a3£2 + ai(£i — T-L)) ± 2^J[T-L — Ci)7Z'^'\ /[a\ + 4£i(?^ — £1)]. Thus, a neces¬ 
sary condition for the trajectory to reach the equatorial plane is that 1-L > £1. 

If the trajectory just touches the equatorial plane but doesn’t go into the lower hemisphere 
then we must have H = L\ and — ^ > £2 > ^. An example of touching the equatorial 
plane is Fig. 21. The different colors in the accompanying graphs correspond to the different 
curves (3.2), (3.3) that make up the trajectories. If the trajectory passes through the equatorial 
plane then we must have % > Ci and 4£f — 4a — 2£i -|- > 0. The angle of crossing lies in the 

interval 

~ “‘^«2£i + a? < cos00 < + ^y^4£2 -4a2£i + af, 

and for each 9q in that interval the possible values of £2 are 


£2 


O3 cos 00 


± \/(£1 — T-L){cos^ 00 £1 — £1 -b ai cos00 


+ 0-2)- 


2 
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Figure 21. Touching the equatorial plane: ai = 3, 02 = 5, 03 


-6, Li = 6, = 6, L 2 = 1. 



Figure 22. Case 1: oi = 3, 02 = 4, 03 = —5, Li = A, L 2 = 0, H = 6. 



Figure 23. Case 1: Oi = 2, 02 = 3, 03 = —4, Li = 3, L2 = 0.5, H = 5. 


Note that for each choice of the constants of the motion, there are always two crossing angles. 
Thus if a trajectory crosses into the lower hemisphere, it must return to the upper hemisphere: 
No trajectory remains confined to the lower hemisphere. Examples are Figs. 22 and 23. If the 
angle is fixed, we can use the formula to get £ 2 - Examples are Figs. 24 and 25. If TZ^ = 0 then 
S{r) = 0. The previous analysis is correct, except the trajectory is a single arc, rather than 
a loop. We call this a metronome orbit. The particle moves back and forth along the arc with 
a fixed period. We give no more details here because we will study the analogous systems in 
our treatment of E16. 
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Figure 24. ai = 8 , 02 = 10, 03 = —10, Li = 40, H = 50, L 2 = —15.1491. 





Figure 25. ai = 8 , 02 = 10, as = —10, Li = 40, H = 50, L 2 = —12.9919. 



Figure 26. ai = 2, 02 = 1, as = —4, Li = 5, H = 8, L 2 = 2. 


Case 2: ai > 02 . Then the portion —1 < si < — 02/01 of the si-axis becomes attractive, the 
trajectories are not closed; each end of the trajectory impacts this interval, assuming TZ'^ > 0 . 
un > Cl then one impact occurs in the northern hemisphere and one impact in the southern 
hemisphere. If 7i < Ci both impacts occur in the northern hemisphere. 

An example for d > Ci is Fig. 26, and for case d < Ci is Fig. 27. If TZ^ = 0 then the 
trajectory is a single arc that impacts the interval — 1 < si < — 02 /ai of the si-axis in the 
northern hemisphere. 
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3.2 System EIQ in Euclidean space as a contraction of S7 

E16 can be described as a Kepler-Coulomb system with barrier in two dimensions [23, 30]. The 
system has Hamiltonian, 


T~L — pi + pi + 


ai 


+ 


022 / 


^ 2 


+ 7/^ x^\/ ^ 


In terms of the rotation generator J = xpy — ypx the defining constants of the motion are PL 
and 


Ci = J^ + 


C 2 = PxJ - 


2 , 022/^2:2 + y 2 y 2 


X^ 

ai y 

2 y^x 2 + 7/2 


+ 03 ^) 


x^ 


- 02 


+ 


c2y^x2 + 7/2 2y^x2 + 7/2 


2 / 

03^. 


The time derivative of a function F of the phase space parameters along a trajectory is given 
dF 
dt 

1 


by F = ^ = {PL, F}, so in Cartesian coordinates we have the equations of motion 


X = 2px 


Px = 


(oix"^ + a27/(3x^ + 2y^) +2a3(x^ + 2/^)^'^^), 


X^(x2 + 7/2)3/2 

1 

y = 2 py, py = -^^—^{aiy-a2). 

In polar coordinates x = r cos 6, y = r sin 6 we have 
P0 = xpy - ypx, rpr = xpx + ypy, 

sin 9 cos 9 

Px = - Pe + cos 9pr, Py = - pe + sm 9pr, 

T 

2 pI oi 02 sin 9 03 ^ 2 sin 9 03 sin2 


PL — Pr - 9 J-1“ 


+ 


J.2 j. J.2 gQg2 Q J.2 (jos2 9 ’ 


— Pe 


+ 


COs2 9 COs2 9 ’ 


, sin 0 \ oi sin 9 02 / sin2 9 1 \ 03 sin 9 

£2 = cos 9pr - Pe Pe -^-1 —vw + 7 ^ I - 


r \ cos 


2 0 2 


r cos2 9 


r ~ j ~ 2 

Setting TZ = {£ 1 , £ 2 } we can verify that 

= 4(£i + — £ 2 ) + o^£i + 2oi02£2 + o^PL. 

Note also that PL = pi + ^ + ■ In polar coordinates the equations of motion are 


(3.4) 




Pr — 2 


£ 1+03 Oi 


J .2 ’ 


r = 2pr, 
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1 


Pe = -3 


cos^ 6 


(02 (cos^ 6 * — 2) — 2a3 sin 0 ). 


From these equations we see that points on a trajectory such that 6 = 0 are characterized by 
Pq = 0 . If such a point has coordinates (6*0, ro,Pro) 0 ) we find that 

^ -02 ± y/al + 4£i(£i + 03 ) ^ 

sin c/Q \ ’ P^o 


ro = 


2(£i + 03) 

02(1+ sin^ 6*0) + 203 sin 60 

(1 — sin^ 6o)(ai sin^o — ^£ 2 ) ’ 


2 -jj -^1 + 03 

Pro=^ - -- 


ro 


'0 


Similarly, if r = 0 then pr = 0 and if such a point has coordinates ( 0 i, ri, we find in 

particular that ri = (oi ± Y/a|~T47Z(£]~+~a3))/(2^). 


Assumption 1 ; oi < 0 , 


02 >0, 03 > 0. 


Possible values taken by the constants of the motion depend on the relative sizes of 02 and 03. 

Case 1 : 03 > 02- Then £1 + 03 > 0 and we must have af + 4 ^(£i + 03) > 0 for trajectories. 
Unbounded trajectories correspond to ^ > 0 , and for "H > 0 perigee occurs at 


Pri = 0 , 

2 £1—02 sin 61 — {Cl + 03) sin^ 9 i 

POi ~ z • 2 a ' 

1 — sm 61 

For 71 = 0 perigee occurs at 

Pri = 0 , 

2 £1—02 sin 61 — {Cl + 03) sin^ 61 

Pei = 


ri = 


Ol + i/ a f + 49 {{Ci + 03) 


sin 01 = — 


ri = 


2 % 

2 ri £2 + 02 
2 £i + 2 o 3 + oiTi 

+ 03 


Ol 


sin 01 = — 


02 


Cl + 03 Ol 


1 - sin^ 01 

Bounded orbits correspond to ^ < 0 , in which case perigee occurs at 

Ol + \/a\ + AT-L{Ci + 03) 


Pri = 0 , 

£1 — 02 sin 01 — (£1 + 03) sin^ 0i 


ri = 


pIi = 


and apogee at 


1 — sin^ 01 


sin 01 = 


2% 

2riC2 + 02 


Pr2 = 0 , 

2 £1—02 sin 02 — (£1 + 03) sin^ 62 

= -r:sin^0^-^ 


r2 = 


2 £i + 203 + oiTi 
Ol - \/ 0,1 + ‘i'H{Ci + 03) 


sin 02 = 


2 % 

2r2C2 + 02 

2 £i + 203 + oir 2 


We must have O2 + 4 £i(£i + 03) > 0 for trajectories. There is an equilibrium point at 


Pro = 0 , peo = 0 , ro = -2 


Ti + 03 

5 

Ol 


sin 00 = 


03 


- \/«i- 


02 


Period of bounded orbit. Since ^ = 2 a/^ _ ^ ^ have 


r dr 


dt 

dt. Thus the time needed to go from perigee ri to apogee r2 is 

fR2 ^ J' 

lim lim / - , = —. 

Rl^ri+R2^r2-J 2 yJ%r‘^ — OiT — (£i + O3) 2 


2^'Hr'^-air-(Ci+ao) 
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Clearly, this is the same time as needed to go from apogee to perigee. Hence the period T is 
T = -^[ lim arctan a — lim arctan B] = -^ ^ . 

2 (_ 7 ^) 3/2 a^-oo 2 (-?^) 3/2 

Thus the period depends only on the energy. 

The trajectories. Eliminating pg, pr in the expressions for T-L, Ci, C 2 we find the orbit 
equations 


A^(sin6») ± 2 Y 5 (sin 6 ') 

r =- . -, N = (0102 - 4a3£2 - 4£i£2)sin0 - 2ai£i - 202 ^ 2 , 

D[sm. 6 ) 

S = (—(jCi + 03) sin^ 6 + Ti — 0,2 sin 9 ^ ( 4 (/ 2 i + a$)(T~l/li — £-2) + o\£\ + 0.2^ T 2aia2>C2) •, 
D = {An{Ci + 03) + O?) sin^ 9 + (4ai£2 + 4a2?^) sin 9 - AUCi + ACl. ( 3 . 5 ) 

Erom structure equation ( 3 . 4 ) and the fact that TZ"^ > 0 we see that the 2 nd factor in S is 
a nonnegative constant. Further, from ( 3 . 4 ) we see that the denominator does not vanish for 
n < 0, that it has a single root sin0 — —2/22/01 for T-i — 0, and a double root for ?/ > 0. 

For the single root case with H = 0 , assuming the plus sign in ( 3 . 5 ), we set z = sin 0 + 2/22/01 
in ( 3 . 5 ) and expand r in a Laurent series series to find the asymptotic expression 

TZ^ (2ai02 — SC1C2 ~ 803/22) 

r = -4^^ H-^- 

ofz^ a\z 

1 + 03) + 16 £|(£i + 03)^ - 8oia2T2(>Ci + 03) + ofa|) ^ 

+ 4 

Note that here TZ^ = — 4 (/ 2 i + az)C\ + a\Ci + 20102/22. If the minus sign is assumed in ( 3 . 5 ), 
then we have the finite limit lim r = — ^^[ 4 / 2 i(/ 2 i + 03) + Oo]. 

Z —^0 

Parametrization of case 1 trajectories. Recall that oi, 02, 03 are fixed constants with 
oi < 0 , 03 > 02 > 0 . To parametrize the trajectories we use perigee coordinates £1, tp, Vp where 
fi = / 2 i + 03 > 0 , — 1 < tp < 1 , Tp > 0 . Then 


LZ + 0 — + —oi. 

rp 

Note that 


.T, ^1 oi tp 1 02 

— “2 H-, £2 — - h - xOitp - -—, 

rp Vp rp I 2rp 

h h 

R > 0 o — > —oi, ?/ = 0 o — = —oi, 


{N + 2 VS) {N - 2 VS) = (aj + 4 (^i - 03)4)T>. ( 3 . 6 ) 

For 7 i < 0 , apogee occurs for radial distance ra = wedge boundaries ^i, S2 

(radial lines containing a point on the trajectory where S vanishes) to exist at all, we must 
have a| + 4 (£i — 03)^1 > 0. Then the righthand boundary will always exist, but the left hand 
boundary will exist only if 2£i > 02 . (However, from the 7Z^ > 0 inequality below, the left hand 
boundary must always exist.) The location of the wedge boundaries is defined by t = sin0, 
where 

_ -02 - + 4^1 (fi - 03) _ -02 + ^/afT 4 £^J£l^^^ 

2 £^ ’ 2 f^ ’ 

(4-a3-02tp--^Hp). 


5i: 
7^2 = 
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ai =-10, a2=3, a3=5, LI =10. L2=-22, H=45 



Figure 28. Case 1 : positive energy. 


In order for the system to correspond to a trajectory, we must have > 0. We will first require 
that 7^ > 0, the generic case. We see that this is violated if either 1) 2^i + oiTp = 0, or 2) tp lies 
outside the open interval (ti,t 2 )- Thus we have the conditions 


y/al + 4£i(4 - 03 ) . 0^ ^^2 + 4^1 (^1 - Q 3 ) 

2£i ^"^2£i 2£i 


2^1 + oirp / 0 . 


For Ti ^ 0 the zeros of the denoiniiicitor occur ut the vulues tof sin 0 such thut 


aiVptp - 02 ± 2 ^{lx + airp)(4 - 03 - a 2 tp - htj) 

2l\ + oirp 

For 77 > 0, if the numerator of the trajectory equation (3.5) vanishes at sin(t) = tn then = t+ 
or t- and ti < tn < t 2 - Since every zero of the denominator is a zero of one of the numerators, 
we must have < t 2 - Further, if 2ii < 02 then D has at most one zero. Thus, if this 

case occurs then there is a single asymptote and a single wedge boundary. 

Summary of necessary Case 1 conditions for trajectories with 7^^ > 0. 


Tp > 0 , 


«2 + Vaj + 4£i(4 - 03 ) 

2h 

For examples, see Figs. 28, 29, 30, 31. 

Case la: TZ^ = 0. In this case equation (3.6) for the trajectories simplifies since 5 = 0. We 
have = (a^ + 4(7i — 03 )^ 1 )^ so 

^ ^ aj + 4(4 - 03)^1 

Because N is linear in t = sin0, it is easily seen that the trajectories are segments of ordinary 
conic sections: lines, ellipses, circles, hyperbolas and parabolas. Now TZ^ = ( ^b+ai^p — 03 — 

02 tp — iitp) = 0, so there are two possibilities: I) 2ii + air^ = 0, and II) — 03 — a 2 tp — = 0. 


1 ) 

2 ) 

3 ) 

4 ) 


— C\ + 03 > 0 , —1 < tp < 1 , 

02 + 4(I’i — a^)(.i > 0, 

2^l + oirp / 0 , 

02 + y/a 2 + 4 £i(£i — 03 ) ^ ^ ^ 

?r?i ^ ^ 
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al = -i, a2=l ,a3=l.l,U=l. L2=l /100, H= -0.1 



Figure 29. Case 1 : negative energy. 



Figure 30. Case 1 : zero energy. 

al a2-l.a3-2, LI -1, L2—l,H-0 



Figure 31. Case 1 : zero energy. 


Case lal): Circular arc. 2.^i + = 0. We can express the constants of the motion in 

terms of tp alone: 


h = 


aiTp 

2 ’ 


n 


2rp 


< 0 , 


C2 = 


Q2 

2rp 
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The equation for the trajectories simplifies to r{9) = 


_ a^+4(£i-a3)£i _ 


N 


= r- 


P) 


so the trajectory is 


a segment of a circle centered at the origin with radius rp. The wedge boundaries are ti = 


- ci 2 — 

airp '■ 


2 ~ airp 


Thus for this case to occur, Vp must satisfy 


O 2 + + 2 aia 3 rp > 0. We conclude that this case always occurs for Vp sufficiently small or 

sufficiently large, but that there is a forbidden intermediate band for which no circular trajectory 
exists. 

Case lall): Conic section arc. — 03 — a 2 tp — £itp = 0. Again we can express the constants 
of the motion in terms of rp, tp alone: 


h 


®3 + 0'2tp 

l-tj 


_ Q 3 + 0-2tp ^ 

r2(l-t2) 


The equation for the trajectories becomes 


T ^2^p) ^2 

rp(l-t2) 2 2rp 


^ ^_ rp{a2 + 2tpa3 + 02 ^ ^ 1 

{2a2tp - ait'^Vp + oiPp + 203) sin^ + (1 - tp)(a2 - aitpVp) asinO + f 3 ' 

where 

{2a2tp - aitlrp + aiVp + 2a3) (1 - t^)(a 2 - aitpPp) 

Q/ = *-^- -■ .. - ^ H = - - -- 

rp{a2 + 2tpa3 + a2tl) ’ rp(a2 + 2tpa3 + a2tp ’ 

and the sign is chosen so that r > 0. Setting y = rsin0, x = rcosO we see that the trajectory 
is a segment of the curve ay + /3r = ±1 or fir = —ay ± 1, so ± 2ay + 1, 

> 0. We have the following possibilities for the curve segments: 


> 0 : ellipse, 0 < : hyperbola, 

/5^ = : parabola, a = 0 : circle, 

/3 = 0: horizontal line. 

The wedge boundaries are tp, —For the circular arc we have rp = —2^^^^^^^. (This 

coincides with Case lal.) For the horizontal line we have tp = < 0, so trajectories are possible 

only for Pp < — 01 / 02 . Only strictly positive energy solutions are possible for the horizontal line. 
For the parabolic arc we have Pp = —203 — 02 — 2 o 2 tp + a 2 tp/oi(l + tp)(l — tp)^. The energy is 
always negative. Only strictly negative energies are possible for the elliptical arc. 

In general, orbits are possible for < 0 if and only if the constants of the motion satisfy the 
inequality 03 + 02 tp/(l — tp)(—oi) < Pp. Trajectories are possible in the positive energy case if 
and only if Pp < 03 + 02 tp/(1 — tp){—ai). 

Case 2: 03 < 02 . Now there are no restrictions on the sign of ii = £1 + 03 . The force in the 
angular direction is always negative. In the first quadrant the force in the y direction is always 
negative. It follows from this that there are no periodic orbits contained entirely in the first 
quadrant. The positive y-axis is repelling, but the negative y-axis is attracting. 

We assume, initially, that TZ^ > 0. Equations (3.5) for the trajectories still hold. The 
condition + 4£i(£i — 03 ) > 0 is now satisfied automatically. There is at most a single wedge 
boundary 


to = sin 00 


-02 + Vaj + 44 {£i - 03) 

2li 


The following general types of behavior occur: 

1. If £1 <0 the force in the radial direction is always negative. The trajectories all impact 
the lower y-axis and are perpendicular to the axis at impact. Since pr is always strictly negative 
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along a trajectory, perigee must occur at yp = —Vp on the negative y-axis. Note that V? is 
linear in T-L, always with positive linear coefficient. The critical value of T-L such that Ti? = 0 
is Eq = [JasTi + JTiTl “ 2aia2T2 — o.\Ci\/[AC‘\ + 02 + The qualitative behavior of the 

trajectories divides into two basic classes, depending on the sign of N{sm8) for 6 = Note 
that 


N{—1) — { 2 C 2 — ai)( 2 £i — 02) + 2ai(o3 — 02), 5 '(— 1 ) — (02 — ■ 

Case 2a: N{—\) > 0. Then if we increase % from Eq to Eq + e for arbitrarily small e > 0 
there will be two points of intersection of the trajectory with the negative y-axis, one on the 
curve r = {N + 2y/S)/D and one on the curve r = {N — 2y/S)/D. The trajectory will remain 
bounded and there will be a wedge boundary. If we further increase "H, leaving £ 1 , £2 unchanged, 
this behavior will persist until the critical value Ei = j[( 2£2 — ai)^]/[(a 2 — 03 )] > Eq, where 
£)(—!) = 0. Here the intersection point with the negative y-axis and the curve r = {N — 2y/S)/D 
has moved to y = —00 whereas the intersection point with the curve r = [N + 2^/S) /D remains 
finite. As T-L is further increased the behavior of the trajectory is that it is unbounded and 
asymptotic to a radial line in the 2 nd quadrant, that there is a wedge boundary and a single 
intersection point with the negative y-axis. 

Case 2b: A"(—1) < 0. Again, we increase T-L from Eq to Eq + e for arbitrarily small e > 0. 
Now there are no points of intersection with the negative y-axis and, in fact, the trajectory 
is non-physical until Ti is increased to e beyond the critical value Ei. The curve is entirely 
described by the function r = {N + 2y/S)/D. The trajectory is unbounded, asymptotic to 
a radial line in the 2nd quadrant and intersects the negative y-axis a single time. There is no 
wedge boundary. Note that for this case, the region {A^(—1) < Ei} is nonphysical. (We 

do not have a complete proof of this but strong numerical evidence.) 

2. If .^1 >0 the behavior is different. The trajectories do not necessarily intersect the negative 
y-axis. Some trajectories are disjoint. The qualitative behavior of the trajectories divides into 
two basic classes, depending on the sign of Y(sin0) for 0 = — |. Note again that 

Y(— 1 ) = {2C2 — oi)( 2 ^i — 02) -I- 201(03 — 02), 'S'(— 1 ) = (02 — aQ)lZ^. 

Case 2c: N{—1) > 0. If we increase % from Eq to Eq + £ for arbitrarily small e > 0 there 
will be two points of intersection of the trajectory with the negative y-axis, one on the curve 
r = [N + 2^/S)/D and one on the curve r = {N — 2y/S)/D. The trajectory will remain bounded 
if < 0 (with a wedge boundary) and be unbounded if > 0. The unbounded trajectory 
will divide into two disconnected parts. If we further increase T-L, leaving £i ,£2 unchanged, 
this behavior will persist until the critical value Ei = \[{2C2 — fli)^][(fl 2 — as)] > Eq, where 
£)(—!) = 0. Here the intersection point with the negative y-axis and the curve r = {N — 2y/S)/D 
has moved to y = —00 whereas the intersection point with the curve r = {N + 2^/S) /D remains 
finite. As Ti is further increased the behavior of the trajectory is that it is unbounded and 
asymptotic to a radial line in the 2 nd quadrant, and a single intersection point with the negative 
y-axis. There is a wedge boundary. 

Case 2d: A^(—1) < 0. Again, we increase Ti from Eq to Eq + e for arbitrarily small e > 0. 
There are no points of intersection with the negative y-axis and, the trajectory is non-physical 
for Ti < 0, but physical for the Ti > 0 interval until T-L is increased to e beyond the critical va¬ 
lue El. Then the trajectory is unbounded and asymptotic to a radial line in the 2nd quadrant. 
Physical trajectories have a wedge boundary. See Fig. 32. 

3. TZ^ = 0. In this case equation (3.6) for the trajectories simplifies since S' = 0. We have 
= (a| -|- 4(^1 — aQ)li)D so 

r(e) - ^ ^2 + 4(^1 - 03)^1 _ ^ al + 4(^1 - 03)^1 _ ^ 1 

N (0102 — 4£2^i)t-|-2(0103 — a2£2 — ai^i) a sin 0 4 -/ 3 ’ 
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al ^-4, a2~3, a3~2, LI »2, L2~-12, H^60 



Figure 32. Case 2d: A^(—1) < 0. 


y = rsin0,x = rcos9 we see that the trajectory is a segment of the curve ay + f3r = ±1 or 
(3r = —ay ± 1, so /3^(x^ + ± 2ay + l,a^ + /3^ > 0. We have the following possibilities 

for the curve segments: 

/3^ > > 0 : ellipse, Q < fP < a^ : hyperbola, 

/5^ = : parabola, a = 0: circle, 

,0 = 0 : horizontal line. 


Since TZ^ = 0 we can express the constants of the motion in terms of a, (3 alone: 


n = - 


a2a^ — aia^ — 2a^a^j3 + a2af3‘^ + ai/3^ 


C2 = - 


2/3 

020 ^ — aia — 2a3a/3 + 02 / 0 ^ 

2 ^ 


ii = 


a2a — ai 


2(3 


Case 2e: circular arc, 0102 — 4 £ 2 ^i = 0. Then T3 = The radius is ro = ±^. There 


- -^24 


ai 


is a wedge boundary at to = sin 0 o = trajectory impacts the negative 

y-axis. See Fig. 33. 

Case 2f: horizontal line segment, 02^2 + = aia^, Vo = ^- The energy is V. = _ 

If the energy is negative there is a wedge boundary at to = sin Oq = ^ 

the trajectory impacts the negative y-axis. If the energy is positive, the trajectory is unbounded 
and there is no wedge boundary, but the trajectory still impacts the negative y-axis. 

Case 2g: parabolic arc. If a = /0 we have 

(02 - asjaf 


C2 = 


01(02 ~ 203 + 2£i) 


n = 


2(2£i-a2) ’ (2£i-a2)2 

and the energy is positive. If a = —/3 we have 


C2 = 


01(02 + 203 — 2£i) 


n = - 


(02 + 03 ) 0 ^ 


2(2^1+02) ’ (2^1+02)2 

and the energy is negative. There is a wedge boundary at to = sin0o = ° 2 +\/a 2 + 4 ^iTi 

the trajectory impacts the negative y-axis. 
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a2~3, a3~2. LI ~-8/5, L2~-15I2, H~-10 



Figure 33. Case 2 e: circular arc. 


Case 2h: elliptic arc. 

_ UiCl - 2 aia 2 C 2 - alih - as) 
al + 4 ii{ii - as) 

the general expression. There is a wedge boundary at to = sin0o = 
trajectory impacts the negative y-axis. 

Case 2i: hyperbolic arc. Here again, 

_ M1C2 — 2aia2C2 — ai(^i — Q3) 

02+ 4£i(^i - as) 

the general expression. There is a wedge boundary at to = sin0o = 
trajectory impacts the negative y-axis. 

There are no radial line trajectories in this case. 


4 The Post—Winternitz superintegrable systems 

These systems are especially interesting as a pair of classical and quantum superintegrable 
systems that do not admit separation of variables, so they cannot be solved by traditional 
methods. We treat only the classical case here and show that its superintegrability allows us to 
determine the trajectories exactly. The quantum case can also be solved. 


4.1 The classical system 

The classical Post-Winternitz system is defined by the Hamiltonian [36] 


1 


^ g (Px + Py) + 


cy 


2 v-^- ' -y/ ' 3 . 2/3 ■ 

The generating constants of the motion are of 3rd and 4th degree: 


£1 = 3plpy + 2pl + 9cx^^^px + 6 


cypy 
3.2/3 ’ 


^ 4 , 4 2 y in 1/3 n 2(9^1^ “ V) 

C 2 =Px + - 12cx / PxPy - 2c -^ 
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Figure 34. c = 5, i? = 2, Li = 1, L 2 = 4. 



Figure 35. c = —12, E = 4, Li = 3, L 2 = 6. 

Here {Li,L 2 } = —188c^, so the constants of the motion generate a Heisenberg algebra. Since 
there is no 2nd degree constant other than Ti, separation of the Hamilton-Jacobi equations is 
not possible in any coordinate system, so separation of variables methods cannot be used to 
compute the trajectories. However, we can find a parametric description of the trajectories. 
Choosing constants % = E, Ci = ii and £2 = -^ 2 , we can verify the two identities 

-p^y + dcx^/'^px -h + QEpy = 0 , 

-i2- ‘^^iPy^ pI - \ {-pI - + QEpyf = 0 . 

From these results we can solve for px and x, y as functions of the parameter py-. 

2{-pi -h + 6Epy) 

Px = , ^=, 

+ 72Ep2 _ 2UiPy + 72E2 _ 13^2 

y = {-‘2^1 - ^QEVy +pI + SIipI - 18Epl + 9 i 2 pl + 72E^ - 18 £ 2 E) , (4.1) 

and equations (4.1) enable us to plot the trajectories. See Figs. 34, 35 and 36. There are no 
closed orbits. The trajectories always impact the y-axis from the left or right, depending on the 
signs of the parameters. 

5 Classical elliptic superintegrable systems 

5.1 Elliptic superintegrability 

Here we extend the ideas in [40, 41] to a system that separates in Jacobi elliptic coordina¬ 
tes [35]. We start with a simple example, the Higgs oscillator S'3 (2.22) [23], expressed in elliptic 
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Figure 36. c = —1, E = 0, Li = 5, L 2 = 5. 


coordinates. 


H = 


1 

sn2(a, k) — sn2(/3, k) 



P} + 


A 

sn2(a, k) 


- - - 1 • 

sn2(/3, k) J 


The Hamilton-Jacobi equation can be taken as H = E for p/s = The separation 

equations in the coordinates a, f3 are 

pl-Esn‘(a,k) + ^^^^+A, = 0, - Esn^(/3, t) + + A, = 0. (5.1) 

Here, H and Li = —p"^ + sv?{a,k)H — are constants of the motion, and Li = Ai on 

a coordinate hypersurface. To find another constant of the motion via the action-angle variable 
construction we introduce M and N which are computed according to 




Pd' 


where Pa,P^ are expressed in terms of coordinates a, /3 by (5.1). Applying action-angle theory, 
see, e.g., [20, 21], we see that M — N is a constant of the motion. The integral M can be written 
as 

sn(Q;, k)da 

^E sn"^(Q:, k) — Ai sn2(Q;, k) — A' 

with a similar integral for N. This can be calculated as an elliptic integral by choosing the new 
variable p = sn^(a, k). We then have the relations 



dp 

da 


2 sn(a, /c)cn(a, /c)dn(a, k), 


cn(a, k) = Y^l — p, dn(a, k) = -\/l — k'^p. 


The change of variables is then tantamount to the computation of 


1 I" cn(Q:, A;)dn(Q;, A:)sn(Q;, A:) da 

2 7 cn(a, A:)dn(a, k)^/E sn‘^{a, k) — Ai sn2(a, k) — A 

^ 1 f __ 

Ak\^ J y/{p - a){p - b){p - c){p- d)' 

where a = 1, 5 = ^ and c, d are roots of ^. A similar calculation applies to the 

substitution p = sn^(/3, k) with exactly the same a, b, c and d. A convenient choice of integral is 


dt 


— a){t — b){t — c){t — d) 
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2 1 / / (fe-d)(t-^ l ib - c){a-d) \ 

y/{a - c){b - d) yY (a - d)(t - 6 ) ’ Y (a - c)(fe - d)y ■ 


This integral can be used for t = p or fi. If we form sn(2A:y^-E(a — c )(6 — d){M — N),k) for 
K = (b — c){a — d)/{a — c)/{b — d) and use the addition formula [35] 


sn(M + u, k) 


sn(u, K)cn(u, K)dn(u, k) + sn(u, K)cn(u, K)dn(?x, k) 
1 — K?SD?‘{u, K)s'n?{v, k) 


we obtain a constant of the motion. In terms of the coordinates it is helpful to write 


Pa 


= VEx 


' {p-c){p-d) 


PP 


= VEx 


' {p-c){p-d) 


Then for Z = we have 

V a—c a—a 


sn(M, K)cn(u, K)dn(u, k) = Z 


p — ab — a 


Vip- c)ip- d) 


p — b p — b 

with a similar expression for sn(u, K)cn{u, K)dn(«, k). We also note that 
2 / x 2 / .2 _ {P - o){P - - c){b - d) 


K^sn(/9, «;)^sn(^, = 


{p — b){p — b){a — c){a — d)^ 


as well as 


y/{p - c){p - d) = 2p^/{l - p){l - k^p)^, 
ViP - c){p - d) = 2py/{l - p){l - k‘^p)^, 


sn(M, k) = 


dn(u, k) = 


'(p-l)(dfc2-l) 

(d-l)(pfc2-l)’ 

/ {p-c){k‘^ - 1 ) 

(c-l)(l-pA:2)’ 


cn(M, k) = 


' (p — d){kP‘ — 1 ) 
(d-l)(l-pA:2)’ 


K = 


'(d-l)(cfc2_l) 

(c — l){dk‘^ — 1 ) ’ 


with similar expressions for sn(v,K), cn{v,K) and dn(r;, k) with p replaced by p. From this we 
deduce that ^^sn{2k-\JE{a — c){b — d){M — N),k) = ^ where JT is a constant of the motion, 

linear in the momenta, viz X = ^ ^\ pPp ~ PP/j,): ™ involution with the 

classical Hamiltonian 


H = -^—{p{l-p){l-k‘^p)pl-p{l-p){l-k^p)pl) - —. 

p — p ^ pp 


This is just the first degree symmetry that we expect for this Hamiltonian. In this case we 
choose the variables p and p to obtain the constant. In particular, cd = —AjE, c + d= —Ai/E. 

This seems to be just a complicated way of deriving a straightforward result. However, we 
can use the TTW trick [40, 41] by looking at the Hamiltonian 


H = 


sn^{pa, k) — sn2(g/3, k) 


pI 


■Pp + 


A 


A 


sn^(pa,k) sn2((jf/3, fc) 
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where p and q are arbitrary nonzero numbers. The separation equations are 


(5.2) 


Here, H and Li = — sv?{pa,k)H + are constants of the motion. To find another 

constant of the motion via the action-angle variable construction we introduce M and N which 
are computed according toM=^f^,N = ^f^, where pa, pp are expressed in terms of 
coordinates a, j3 by (5.2). Again M — N is a constant of the motion. The integral M can be 
written as 

If sn(»a, k) da 

M = - , ^ J 

2 J E sn‘^{pa, k) — Ai sn‘^{pa, k) — A 

with a similar integral for N. 

This can be calculated as an elliptic integral by choosing the new variable p = ST?{pa,k). 
We then have the relations 


= 2p sn(pa, A:)cn(pa, /c)dn(pa, fc), 
da 

cn{pa, k) = sjx — p^ dn(pa, k) = \/l — k'^p. 


The change of variables is then tantamount to the computation of 

cn(d, /c)dn(d, A;)sn(Q;, k) da 
cn(Q;, A:)dn(Q!, 

^ 1 f __ 

Ak\^p J {p - a){p - b){p - c){p- d)' 

where a = pa^ a = h = ^ and c,d are roots oi p^ — ^p — A similar calculation applies 
to the substitution p = sv?{qj3,k) with exactly the same a, b, c and d. A convenient choice of 
integral is 



dt 


y/ {t — a){t — b){t — c){t — d) 


y/{a- c){b - d) 


sn 


-1 


{b — d){t — a) (b — c) (a — d) 

{a — d){t — b)^ {a — c){b — d) 


This integral can be used for t = p oi p. If we now form sn(2A:Y^£'(a — c)(6 — d)pq{M — N)) 
and use the addition formula 


sn(gu + pv, k) 


sn{qu, «))cn(pu, K)dn{pv, k) + sn{pv, K)cn{qu, K)dn{qu, k) 
1 — K‘^sn^{qu, K)sn^{pv, n) 


we obtain a constant of the motion. Thus computation yields a rational constant of the motion 
whenever p/q is a, rational number. 

Let us restrict ourselves to p = 1, q = 2. Then u = ajl and v = f3 where 


Pa 


= VEy 


' {p-c){p-d) 


PP 


= VE\ 


{p-c){p-d) 

P 


Y = 


b-d 
a — d 
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We form 

Ti = sn( 2 rt, K)cn(ti, K)dn(ti, k) 

^ _ 2 y _ \/(P -a)ip- b)ip -c){p- d){^i - c)(^ - d) _ 

{fi — b){p^{a + h — c — d) + 2p{cd — ah) — cd{a + 6 ) + ab{c + d)) ’ 

T 2 = sn(u, K)cn( 2 ti, K)dn( 2 r;, k) 

= —YA ——- ((c + 6 — a — d)p‘^ + 2 {ad — bc)p + cb{a + d) — ad{b + c)) 

\ p — 0 

X ((d + b — a — c)p^ + 2{ac — bd)p + cb{d — a) + ad{b — c)) 

X ((a + b — c — d)p^ + 2{cd — ah)p + ah{c + d) — cd{a + 6 )) ^ . 

We also note that 

T 3 = K^sn( 2 rt, K)‘^sn{p, 

^ 4 _ id -a)ip- a){p -b){p- c){p -d){b- d){b - c) _ 

{p — b) ((a + b — c — d)p‘^ + 2{cd — ab)p + ab{c + d) — cd{a + 6 ))^ ’ 

ViP - c){p - d) = 4^V(l-/i)(l-/cV)^, 

V(P - c){p - d) = 2 py'{1 - p){l - k^p)^, 

and p = sn^(a, k) and p = sn^(2/3, k). From this we deduce that ^sn(4/cY^E(a — c){b — d){M — 
N), k) = P, where P is a constant of the motion rational in the momenta. We can calculate the 
components of this constant as follows: Ti = ni/di, where 

m = -16(1 - p)[l- k'^p)ppk'‘^y/l - pppPf,, 

di = {k^p - l){{k^ki - P(1 + k2))p^ + 2{Ak^ + E)p - Ai - A{k2 + 1)), 

T 2 = n 2 /d 2 , where 

n2 = aJI - p[-{E‘^k'^ + k^{Kl + AAE))p^ + {-k^{ - 8AE + 2AiP - 2A?) 

+ k^{-8AE - 4AiP - A?) + 2AiP)p=^ + 6{2k^AE + AE + k^{kl + AE))p^ 

+ {-2k‘^kiA + 2k^{-kl + 2kiA - AAE) - 2A? - 2AiA - 8 AE)p 
+ {-k*A2 + AAE + {2k^ - 1)^2 - A?)], 
d2 = Vk^P-^[{k^^i - P(1 + k^))p^ + 2{Ak^ + E)p - Ai - A{k^ + l)]^ 

Ts = ns/da, where 

ns = 4(p - 1) (ypk^ -l){p- 1) [-Ep^ + kip + A) {-E + + AiA:^), 

ds = {k‘^p - 1) [(A^Ai - P(l + k^))p^ + 2{Ak^ + E)p-ki- A{k^ + l)]l 

We could, of course, use the expressions above given in terms of c + d and cd. 

If we set X = sn^(a, k) = p and y = sn^(2/3, k) = p then H and Li become 

H= -^^{x{l-x){l-k‘^x)pl-iy{l-y){l-k‘^y)pl) - —, 

X — y^ \ / V ' y' xy 

(1 - k^x)pl - 4(1 - y) (1 - yk^)pl) - A 0 + 


(5.3) 



42 


Y. Chen, E.G. Kalnins, Q. Li and W. Miller Jr. 


For ^ = 0 the extra constant of the motion is G = N/D where 

N = y/y{{-yxk^ - x^k -x"^ + 2x- y)pl - 4y{y - 1)(-1 + k‘^y)pl 
+ 4x(y -!)(-! + k^y)pxPy) , 

D = \/y - l[[x^k^ + x'^yk^ - x^ - 2xyk^ + y)pl + 4y{y - 1)(-1 + k^y)py 

- 4y(-l + k^y){x - l)pxPy)- 

Setting R = {Li, G} we obtain the polynomial structure relations 

- 4(1 - k'^G^) (Li + G^H - LiG^) = 0, 

{G,ii} = 2{G^ - l){G^k'^ - 1 ), 

{Li, ii} = -4G((2 - k‘^)Li + 2k'^HG^ - 2k'‘^LiG‘^ - H). 

Turning to the case when M 7 ^ 0 we have G = N/D, where 

N = ^,/y — 1 [x‘^y[x‘^yk‘^ + — x^ — 2xyk^ + y)p\ + 4x^y^(y — 1)(—1 + k‘^y)py 

- 4x^y^(-l + k^y)pxPy + A{x - yf ], 

D = \/—1 + k‘^y[x‘^y[—2xy + y — k^x'^ + x^yk^ + x^)p\ + 4x^y^(y — 1)(—1 + k‘^y)py 

- Ax^y^{y - 1) (-1 + k^x)pxPy + A{x - yf]. 

The structure relations are 

{G,R} = -2 (g2 - l){k'‘^G‘^ - 1 ), 

{Li, ii} = -AG[2Lik^G‘^ + 2Ak^G‘^ - 2HG^ - h - hk^ - 2Ak^ + 2H], (5.4) 

R^ + A[Ak*G‘^ - Lik'^G^ + LiG^ - 2Ak‘^G + A - LiG^ + 2HG‘^ + Li-H -H^] = 0 . 

All of the examples that we have constructed using the addition theorem for elliptic functions 
obey polynomial structure equations, but we have no proof that this is true in general. 


5.2 Example 


H = 
Li = 


A 


— X 


l((l - ^")Px - 4(1 - y^)pl) - 

D2 A{x- + - 2) 

(y2 _ a.2) ^ (1 _ 3.2)(1 _ y2) ’ 


L 2 



2x^y + y^)pl + {-Ax'^y + Ax^y^)pl + (-dx^y^ + Ax^)pxPy - 

- 2y2x2 + y2)p2 (4a™2 _ 4a;2y2)p2 (_4a;y + Axy^)pxPy + 


This is the special case of (5.3), (5.4) where we have replaced x, y by 1 — x^, 1 — y^, respectively, 
set A: = 0 and rescaled. Note that H = 4y^ + + (T^^’ ~ ^Py) ^ 

- (f-x2)'(Ly2)2 ■ We have {H, Li] = {H, L 2 ] = 0, and with R = {Li, L2}, the polynomial 
structure relations are 


{Li, i?} - I 6 L 2 T 1 + 2>2HL2 - 2>2LIH = 0 , {L2,R] - 8 (l - L|) = 0 , 

R^ + I6L1 - IQLILi + 32HLI - IQL^H - 16H + 16A = 0. 
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Using the three constants of the motion to eliminate the momenta terms we find 
2 _H{x^ - 1)^ + Li(x^ -1) - A 2 _ “ 1)^ + ^liy^ - 1) - ^ 

- (l-x 2)2 ’ Py ~ 4 (l-y 2)2 ’ 

and the trajectories must satisfy the implicit equation 

0 = (^Hx'^iA + Li - Lix^ - Hx^ + 2 Hx‘^ - H)y^ + ^ {-H + Li + Hx^ + A)^^ L| 

+ (sHx'^Li - \h‘^x^ - + aH- 3 H^x^ + Lix'^A + 2H^x^ + ^Hx^A 

- ALi - LiHx^ + LiH + 2 H‘^x‘^ - -L\ + - 3 Lix‘^H - 2 Hx^a\ 

2 2 / 


vL2 


+ ^{-H + Li + Hx'^ + Afy'^ 

-x‘^{A + Li- H) {A + Li- Lix^ - Hx^ + 2Hx‘^ - H). 


The left-hand side of this equation is quadratic in y, so we can solve for y as a function of x. 
We obtain the solutions y = where 

N = L 2(-1 + xf{x + - 2L2{-2 Ax‘^ - 3Lix^ - Lix® + 3Lix^ + 3Ax^ + A + Li)H 

-j- L2(A -|- Li) [a -|- L\ — 2Lix‘^'), 

d2 

S = -xHA + Li - Lix^ - Hx^ + 2Hx'^ - H) i-H + Li + Hx^ + A)—, 

\ / \ ^16 

D = —(—1 -I- x)‘^{x + l)^(x^ -|- 2L2X + 1 ) {—x^ -|- 2L2X — 1 )H^ 

+ {ALlAx'^ -h 2Lia;^ - 2Li -b dL^Lix^ -2A- 4L^Lix^ -b 2Ax'^)H + {A + Lif. 


Let (a;o,yo) be a point on the trajectory such that = 0. A straightforward calculation gives 


2 2H -Li± JL'i + AAH 
xo = - ^ -> yo = -L2, 

2 A- H - HLi -b 2HL^ + Li - L^Ll 
4(L2-1)2(L2 + 1)2 

Assumption: A > 0. Note that 


/ -Li± v/L2+4Ag \ 

V )' 


(5.5) 


It follows that if we are in the region x > 1 and H > 0, then since y^L^ -b AAH > |Li|, only 
the plus sign is possible in (5.5). We conclude that the trajectory must be unbounded to the 
right. Examples are Figs. 37 and 38. (The different colors in the figures correspond to different 
functions describing the trajectories.) 

We investigate what happens if the trajectory intersects the lines y = Lx. If the intersection 
occurs for y = x then we find that x = —L 2 or px = 2py, or both. If the intersection occurs for 
y = —X then we find that x = L 2 , or px = —2py, or both. In both cases, we have = 4p2 = 

• Thus, in both cases if IL 2 I 7 ^ 1, the trajectory intersects the lines y = Lx 

at finite momentum. However, the velocity is infinite. Indeed x = {H, x} = 2px{l—x‘^)/{y‘^—x‘^), 
which blows up unless p^ = 0. This suggests that there are no bounded orbits in the square 
— 1 < X, y < 1. However we find some bounded orbits in the region x > 1 for negative energies. 
Examples of these bounded orbits are Figs. 39 and 40. 
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Figure 37. Case: A = 0.5, Li = 0.1, La = -0.4, H = 1. 



Figure 38. Case: A = 1, Li = 0.2, L 2 = 0, H — 2. 



Figure 39. Case: A = 5, Li = 10, La = 3, 74 = —1. 



Figure 40. Case: A = 1, Li = 10, La = —2, H = —2. 


5.3 Another version of the example 

We set y = iz where z is real, multiply H hy —1 and L2 by —i to get the system 

A 


(1 4(1 + ^2)^2 

ri — - 5 -h 


(1 — x^)(l + z‘^) ’ 


+ z'^ 
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{I + z^){l — + Apl) A{x^ — z^ — 2) 

(1 + z 2 )(l — x^) ’ 


L 2 


4 9 

zx^p^ 


Ax^PzPxZ^ — ‘^x^PzPx — 2x‘^zp‘^ + ix'^plz + 4 x ^ p ^2;^ — p'^z^ - 
4p2x2 — ApxZPzX — ^p‘^x‘^z‘^ + 2p%z‘^x‘^ — p^z"^ — 4:XZ^PzPx + 


Ax^z(z^+x^)^ 
■ (1+22)(1-x2)2 

Aa; 2 (a; 2 + 22)2 

(1+22)(1-x2)2) 


We will consider this system for ^ > 0 and resticted to the vertical strip of points (x, z) where 
— 1 < X < 1. Then H must be > 0. The motion is unbounded, there are no bounded orbits, and 
the origin ( 0 , 0 ) and the boundary lines x = ±1 are repulsive. 

The polynomial structure relations are R = {Li,L 2 } and 


{L 2 , 4?} + 8^1 + T 2 ) — 0) {-bi, — 16T2.bi — 32 HL 2 — 32 L 2 H — 0, 

-R'^ + 16A + I 6 L 1 + 16L|Li + 32HLI + IGL^H + 16H = 0. 


Expressing the momenta in terms of the coordinates we have 

2 _ {A + H + x^H + Li- x^Li - 2Hx^ 

Px - (l-x2)2 ’ 

2 _ 1 {2Hz‘^ + Hz^ + Liz2 + A + H + L 1 ) 

- 4 0+^^? ■ 

The equations for the trajectories turn out to be y = , where 

N = L2{1- x^Yh"^ + 2L2(Li + 3^x^ + A - Lix® + 3Lix^ - 2Ax^ - 3 LixYH 

+ L2 (Ti + A) (Li -j- a — 2Lix^), 

2 

S = - — (A + H + Hx^ + Li- Lix^ - 2Hx^) (A + H + Li- Hx^fR^, 

16 

D = (1 - x‘^f{x^ + 4.Llx^ + 2 x^ + 

+ (-2Ax^ + 2Li - 4L|Lix^ + 2A - 2Lix^ + 4L|Lix2 + ALIAx^H + (Li + Af. 


The points {xq, zo,pxQ,Pzo) on the trajectory where pxQ = 0 satisfy 

2 f + i?(l + + Li + L1L2) 

Zo = -L 2 , Pzo = l (1 + L2)2 ’ 

2 El ± VE2 - AAH 
^ 0 - 1 + • 

Requiring first that R? > 0, we see that Li < 0 is a necessary condition for trajectories (otherwise 
5 < 0 in the strip — 1 < x < 1. If \Li\/2H > 1 then a further necessary condition for S > 0 
in a strip is Li H A A < 0. If |Ei|/2fI < 1 then a necessary condition is L2 — AAH > 0. 
Examples here are Eigs. 41 and 42. Note that here there are 4 different trajectories corresponding 
to a single choice of constants of the motion. 

Case when R^ = 0. When Li = —9, L 2 = 0, ^ = 1, = 8 , we have R^ = 0 such that S = 0. 

The trajectory in this case is a straight line on the x-axis. The two boundary points where Pa; = 0 
for X > 0 are 0 and 0.9354. The velocity of the particle is x = {H,x} = ^ • Substituting 

z = 0 and the expression for px in terms of x into the equation gives,x = \J~ ^ Thus 

lim X = 00 , so the particle will go through the origin instead of bouncing back. 
x^O 
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Figure 41. Case: A = 1, Li = —5, L 2 = 3, = 2. 



Figure 42. Case: A = 1, Li = —10, L 2 = —3, H = 10. 


A Review of some basic concepts in Hamiltonian mechanics 

Hamiltonian mechanics is a reformulation of classical Newtonian mechanics alternate to La- 
grangian mechanics. In the Hamiltonian formalism a physical system describing the motion of 
a particle at time t involves n generalized coordinates qj(t), and n generalized momenta Pj{t). 
The phase space of the system is described by points {pj,qj) G The generalized coordina¬ 
tes q'jS and generalized momenta p' s are derived from the Lagrangian formulation. The dynamics 
of the system are given by Hamilton’s equations [1, 12] 

dqi^ dn dpi^_dn 

dt dpj' dt %’ ^ ’ 

where % = T-L{qi,... ,qn,Pi, ■ ■ ■ ,Pn,t) is the Hamiltonian, which in this paper corresponds to 
the total (time-independent) energy of the system: T-L = T + V , where T and V are kinetic and 
potential energy, respectively. Solutions of these equations give the trajectories of the system. 
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Here T is a function of q alone while E is a function of q alone. Explicitly, 

^ + ^(q), 

j,k 


(A.2) 


where is a contravariant metric tensor on some real or complex Riemannian manifold. That 

n 

is g~^ = det{g^^) / 0, g^’^ = g^^ and the metric on the manifold is given by = Y1 9 jkdq^dq’^, 

j,k=l 

where (gjk) is the covariant metric tensor, the matrix inverse to (g^^)- Under a local transfor¬ 
mation g'- = /; (q) the contravariant tensor and momenta transform according to 


i9T = }^7^7^ 


j,k 


jk 

dqj dqk ^ 


dqj 


n 

j=l t 


SO % is coordinate independent. Here m is a scaling parameter that can often be interpreted 
as the mass of the particle. Eor the Hamiltonian (A.2) the relation between momenta and 

n n 

the velocities is pj = mYl djiqe, so that T = g^'^{^)PjPk = t E 5a(q)%4- Once the 

i=l j,k e,h=l 

velocities are given, the momenta are scaled linearly in m. In mechanics the exact value of m may 
be important, but for our purposes can be scaled to any nonzero value using the above formulas. 
To make direct contact with mechanics we will often set m = 1; for structure calculations we 
will usually set m = 1/2. 

The Poisson bracket of two arbitrary functions ^(p, q), H(p, q) on the phase space is the 
function 



dA dB 

dpj dqj 


dA d]3\ 

dqj dpj ) 


The Poisson bracket obeys the following properties, for A, B, C functions on the phase space 
and a, b constants. 


anti-symmetry: {A,B} = —{B,A}, 

Bilinearity: {^, mB + nC] = m{A, B} + n{A, C}, 

Jacobi identity: {^, {B, C}} -|- {B, {C, ^}} -|- {C, {A, B}} = 0, 
Leibniz rule: {A,BC} = {A,B}C + B{A,C'\, 
chain rule: {f{A),B} = f'{A){A,B}. 


In terms of the Kronecker delta 6jk, coordinates (q, p) satisfy canonical relations {pj,Pk} = 
{qj^qk} = 0, and {pj,qk} = djfc- Using the Poisson bracket, we can rewrite Hamilton’s equa¬ 
tions (A.l) as 


{'H^qj} 


dt 


dn 

dpj' 


{'H,Pj} 


d^ 

dt 


ddi 

dqj' 


Eor any function T'(q, p), its dynamics along a trajectory q(t), p(t) is Thus 

will be constant along a trajectory if and only = 0. If {TijiF} = 0, we say 

that T" is a constant of the motion. 


A.l Classical integrability 

A system with Hamiltonian PL is integrable if it admits n constants of the motion Ci = PL, 
£ 2 ) • • •) that are in involution: 

{£j,£fc} = 0, l<j,k<n, 


(A.3) 
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and are functionally independent in the sense that det (§^) / 0. Suppose % is integrable 
with associated constants of the motion Cj. Then by the inverse function theorem we can 
solve the n equations £j(q, p) = Cj for the momenta to obtain = pfc(q, c), k = 
where c = (ci,...,Cn,) is a vector of constants. For an integrable system, if a particle with 
position q lies on the common intersection of the hypersurfaces Cj = Cj for constants Cj, then its 
momentum p is completely determined. Also, if a particle following a trajectory of an integrable 
system lies on the common intersection of the hypersurfaces Cj = Cj at time to, where the Cj 
are constants of the motion, then it lies on the same common intersection for all t near to- 
Considering Pj(q, c) and using the conditions (A.3) and the chain rule it is straightforward to 
verify ^ Therefore, there exists a function u{q, c) such that p£ = ^, £ = 1,... ,n. Note 

that Vjiq, = Cj, j = 1,... ,n, and, in particular, u satisfies the Hamilton-Jacobi equation 

where E = ci. By construction det( g^^g^^ ) / 0, and such a solution of the Hamilton-Jacobi 
equation depending nontrivially on n parameters c is called a complete integral. This argument 
is reversible: a complete integral of (A.4) determines n constants of the motion in involution, 
Vi,..., Vn- 

Theorem A.l. A system is integrable if and only if (A.4) admits a complete integral. 

A powerful method for demonstrating that a system is integrable is to exhibit a complete 
integral by using additive separation of variables. 

It is a standard result in classical mechanics that for an integrable system one can integrate 
Hamilton’s equations and obtain the trajectories [1, 12]. The Hamiltonian formalism is is well 
suited to exploiting symmetries of the system and an important tool in laying the framework 
for quantum mechanics. 

A.2 Classical superintegrability 

Let = (/i(q, p),... ,/ 7 v(q)P)) be a set of N functions defined and locally analytic in some 
region of a 2n-dimensional phase space. We say E is functionally independent if the N x 2n 
matrix rank N throughout the region. Necessarily N < 2n. The set is functionally 

dependent if the rank is strictly less than N on the region. In this case there is a nonzero analytic 
function T of A variables such that F{fi,..., fn) = 0 identically on the region. Conversely, 
if F exists then the rank of the matrix is < N. We say that the Hamiltonian system H is 
(polynomially) integrable if there exists a set of A = n constants of the motion Ci =%,... ,Cn, 
each polynomial in the momenta globally defined (except possibly for isolated singularities), 
that is functionally independent and in involution: {Cj,Ck}=£) for 1 < j, A; < n. A classical 
Hamiltonian system in n dimensions is maximally {polynomially) superintegrable if it admits 
2 n — 1 functionally independent, globally defined constants, polynomial in the momenta (the 
maximum number possible [33]). At most n functionally independent constants of the motion 
can be in mutual involution [I]. However, several distinct n-subsets of the 2n — 1 polynomial 
constants of the motion for a superintegrable system could be in involution. In that case the 
system is multi-integrable. An important feature of superintegrable systems is that the orbits 
traced out by the trajectories can be determined algebraically, without the need for integration. 
Along any trajectory each of the symmetries is constant: Cg = Cg, s = 1, ..., 2n — 1. Each 
equation /^^(q, p) = Cg determines a (2n — l)-dimensional hypersurface in the 2n-dimensional 
phase space, and the trajectory must lie in that hypersurface. Thus the trajectory lies in the 
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common intersection of 2n — 1 independent hypersurfaces; hence it must be a curve. Another 
important feature of the trajectories is that all bounded orbits are periodic [34], 

The polynomial constants of the motion for a system with Hamiltonian % form the (poly¬ 
nomial) symmetry algebra 5-^ of the system, closed under scalar multiplication and addition, 
multiplication and the Poisson bracket: Indeed if H is a Hamiltonian with constants of the 
motion £, fC. Then aC + /3/C, CK. and {T,/C} are also constants of the motion. The n defining 
constants of the motion of a polynomially integrable system do not generate a very interesting 
symmetry algebra, because all Poisson brackets vanish. However, for the 2n — 1 generators of 
a polynomial superintegrable system the brackets cannot all vanish and the symmetry algebra 
has nontrivial structure. The degree (or order) 0{C) of a polynomial constant of the motion L 
is its degree as a polynomial in the momenta. Here P/ has degree 2. The degree 0{Fk) of a set 
of generators = {£i,...,£fc}) is the maximum degree of the generators. Let S = Sp^. be 
a symmetry algebra of a Hamiltonian system generated by the set F^- Clearly, many different 
sets can generate the same symmetry algebra. Among all these there will be a set of genera¬ 
tors Fj^^ for which i = 0{F^^) is a minimum. Here £ is unique, although Fj^^ is not. We define 
the degree (or order) of S to be £. The theory of 2nd degree superintegrable systems has been 
developed in papers such as [8, 16, 17, 18, 19, 22, 28, 37]; all such systems are known as are the 
structures of the symmetry algebras, all of which close at degree 6. 
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